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CHANNEL-SLOPE FACTOR FLOOD-FREQUENCY ANALYSIS 


Manuel A.M. ASCE 


ABSTRACT 


Annual flood peaks New England have been related many hydrologic 
factors. The main-channel slope has been found next importance 
drainage-area size. The slope for that part the main channel between 
and percent the total distance above the gaging point provides the best 
correlation with flood magnitudes. 


Many agencies and individuals are studying the causes floods and their 
relation known meteorologic and physiographic conditions, and are develop- 
ing generalized flood-frequency relations for the purposes hydraulic de- 
sign, flood insurance, flood-plain occupancy, and other uses. 

Data for New England have been used the first phase such study 
the Geological Survey. The records runoff and precipitation for this 
region are numerous and long records for any other region the 
United States; topographic quadrangle maps covering the entire area are 
available; historical flood data are available, dating back the earliest 
times settlement. 

Peak discharge records 170 gaging stations were used the analysis. 
These records represented all the stations New England having least 
years record, which the effect regulation annual peak discharges 
was known, judged, negligible. These stations were located through- 
out the region, though not densely the northern part. The range the 
size drainage areas was from 1.64 9,661 square miles. total 
170 station, stations drained areas less than square miles; stations 
drained areas less than square miles; stations drained areas less than 
100 square miles; and 106 stations drained areas less than 200 square miles. 

practice, there are many ways converting set flood data 
single station flood-frequency curve relation for that station, 
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determine the floods various recurrence intervals. was considered that 
the use mathematically fitted curves involving prior assumptions 
distribution might obscure the very relationships that were being sought, and 
that graphically drawn curves, conforming closely possible the 
original data each station, would most desirable. The plotting position 


each flood record was computed using the formula 


the number years record and the rank the flood, starting with 
the highest. This formula has theoretical basis and also agrees closely 
with the length the period record for the top flood, thus again conform- 
ing closely the basic data. 

From the frequency curves, values were selected for the 1.2-, 2-, 5-, 10-, 
and 25-year floods 170 stations, and for the 50-year flood stations. 
Where the length record was less than years, correlation with nearby 
stations was used ascertain the appropriate plotting positions within 25- 
year period for the annual peaks record. The floods for each the stated 
recurrence intervals were then separately correlated with large number 
physiographic and meteorologic factors. Initially, the relations were ana- 
lyzed graphical multiple-correlation techniques combined with rank- 
correlation methods. The following parameters were studied: 


Drainage area 


Slope 
Channel slope No. 
Channel slope No. 
Average land slope 
Tributary channel slope 


Water-supply index 
Mean annual runoff 
Mean annual precipitation 
Maximum 24-hour precipitation 
10-year 60-minute intensity 
Total rainfall volume large storms 


Shape factor 
L2/A, equivalent basin length divided basin width 
L/A, basin length divided drainage area 
the summation small subdivisions the drainage area, each 
multiplied the distance travel the gaging point. 


Mean elevation 

Percentage lakes, swamps and reservoirs 
Stream density 

Temperature (January average) 


Most the physical characteristics were computed described 
Langbein [1947]. “Channel slope No. one such characteristic. the 
mean slope all parts channel draining least percent the total 
drainage area,—thus may include the principal tributary streams. 
“Channel slope No. the same that described Bigwood and Thomas 
[1955]; the weighted main-channel slope averaged with the tributary 
channel slopes described these authors. 
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The mean annual precipitation for each drainage area was supplied the 
Weather Bureau computed the preparation Hydrologic Atlas 
and Nordensen]. Other precipitation factors were based vari- 
ous Weather Bureau publications Geological Survey flood reports. 

Many the hydrologic factors, such “channel slope No. “stream 
density” were available from Water Supply Paper 968-C, Langbein, [1947], 
though not for all the 170 stations usable New England. Because the 
labor involved computing these, was considered satisfactory, the 
preliminary phases the study, compare the various factors and test 
their significance using only the stations for which most the 
hydrologic factors were already evaluated. 

the graphical procedure used, Qn, the n-year flood, was first plotted 
(on log-log paper) against the drainage area. Through the plotted points, 
average curve was drawn, and the residuals from this curve were determined 
for each station and ranked order magnitude. The values for all other 
parameters were then ranked and the differences between ranks were used 
compute the Spearman rank correlation coefficients[ Wallis and Roberts, 
1956]. The parameter with the highest correlation (considering only one 
from each numbered group the preceding table) was then used the next 
variable the graphical correlation. After the proper adjustments, plot 
(adjusted) first residuals versus the first parameter values led set 
secondary residuals, which was turn tested against the remaining parame- 
ters. 

Table shows the results graphical correlation for the 2-, 10-, and 
25-year floods for all significant parameters. This shows the decrease 
the standard error when successive independent variables were included. 
each case, after independent variables were used, further variables 
did not appreciably decrease the standard error. notable that factors 
evaluating tributary slope alone were not found significant, once the 
effect the main channel slope, expressed channel-slope factor no. 
was accounted for. The graphical standard errors were determined the 
range within which two-thirds the plotted points lie. Actually these are 
not standard errors estimate, but are approximations adequate for these 
exploratory studies. The standard errors are shown terms the original 
log units and the equivalent percentages these represent. One log unit 
equals one log cycle and decimal part log unit that same 
part the linear distance representing one log cycle. 

would expected, the standard error increases with the increase 
the size flood. will noted that for each size flood, the standard 
error with drainage area (A) drops sharply after channel slope (S) in- 
cluded, and improved only small amount thereafter. The average 
percentage improvement with 22, 38, and percent respectively, 
whereas the other parameters after and show combined improvements 
not exceeding percent. 

will noted also that water supply indices, show but small 
effect. The relatively minor influence rainfall attributed the uniform 
pattern precipitation throughout New England and the unsatisfactory 
rainfall parameters far available. 

The importance channel slope shown the large percentage the 
error which accounts for. However there has been unique universal- 
accepted way evaluating channel slope. Some hydrologists have used the 
total drop from the head the longest watercourse. Others have used 
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Table 1.--Results Graphical Multiple Correlation 


Standard error 


log percent 
units ave 


A=Drainage area; S=Channel slope No. R=Mean annual runoff; 
P=Mean annual precipitation; 60-minute maximum 
precipitation; L/A=Shape factor. 


weighting methods which evaluate the slope all along the main channel. Some 
have used parts the main channel, such the lower three-quarters. Still 
others have some manner combined the slopes tributary streams with 
the main channel. Part the difficulty that methods defining the main 
channel differentiating the tributary streams are not fixed. 

Channel slope No. was found the best slope factor those evalu- 
ated. Its computation involves the determination the points the larger 
streams which percent the total area drained. This laborious 
task, particularly regional study where may have done for 
hundreds stations. was therefore considered desirable attempt find 
simpler, and possibly better channel slope factor. was also considered 
desirable separate the effects tributary streams from those the main 
stream. was decided (1) use only the main channel channel-slope 
parameter, leaving the tributary slopes for separate consideration, (2) de- 
fine the main channel, above each stream junction, that channel draining 
the largest area, (3) make exhaustive study find way expressing 
the main-channel slope that most closely related peak discharge. 

Topographic maps were used develop channel profiles for all 170 sta- 
tions New England. the upstream end each profile was extended the 
drainage divide beyond the end the stream shown the topographic map. 
Distances were measured from the gaging point each contour crossing (ex- 
cept where these were too dense) and channel profiles were plotted from 
these figures. 
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was considered that the most upstream part the stream, the steep 
headwaters, might affect the slope out proportion the volume water 
furnished the headwater area. the downstream end the slope might not 
indicative that affecting the size peak discharges because gaging 
sites are usually located riffles natural controls. was therefore 
postulated that the part the main channel whose slope would best corre- 
late with peak discharge would the one excluding the extreme headwater 
reach, and possibly some the extreme downstream end. 

The total lengths the streams were divided into tenths below 0.7 the 
distance from the gage and half-tenths beyond that. The object was com- 
pute slopes for all possible portions the channel, for example, 0.95 0.0, 
0.7 0.1, 0.6 0.4, etc., the figures representing each case the fraction 
the total distance upstream from the gage. There were such combina- 
tions. Two additional slope factors were computed. These were the con- 
stants the regression equation relating the logarithm the rise the 
logarithm the distance from the gage, equivalent drawing straight line 
through the channel profile plotted logarithmic scale. These two factors 
represent integrated slope for the entire main channel, rather than mere- 
the slope between two points. 

determine the optimum slope factor, multiple correlation analyses 
were made with the 1.2-, 2.33-, 5-, 10-, 25-, and 50-year floods de- 
pendent variables and drainage area and each the slope factors the 
independent variables. The change from the 2-year the 2.33-year flood 
was made this time order use value corresponding the mean an- 
nual flood common use; results are comparable. The best slope factor 
would defined the point minimum standard error. The correlation 
analyses were such prodigious task that the data were prepared for auto- 
matic solution digital computer, the Datatron 205. The programming 
covered computation all the slopes, the logarithmic slope factors, the 
standard errors with drainage areas alone, and the standard errors and cor- 
relation coefficients with drainage and each the slope factors. Data for 
170 stations was used these correlation studies. 

Results the Datatron solutions are shown Figures and These 
figures show contours representing equal values the standard error. 
Figure for example, the standard error using the slope between points 0.5 
and 0.4 the total distance above the gage 0.193 log units. Figure shows 
all the standard errors which the contours are based whereas Figure 
shows only the contours. Note the similar patterns shown the controus 
for each size flood and the minimum the upper left each figure. 
weighting process reveals that the slope between points and percent 
above the gage would satisfactorily give the minimum standard error for all 
floods, with little accuracy lost for any particular size flood. Standard 
errors using the two logarithmic slope factors were higher for all floods 
than those using the percent slope. 

Table shows the standard errors (in log units and average percent) for 
each size flood, using drainage area and the “85-10” percent slope in- 
dependent variables. The improvement each can seen when slope in- 
troduced. regularity shows the percentage the error explained 
slope, which close percent each case. The standard error for the 
50-year flood seems inconsistently small and the contour pattern some- 
what erratic. The irregularities may rise because the 50-year flood based 
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UPSTREAM POINT 


DOWNSTREAM POINT PROFILE 


FIG. VARIATION STANDARD ERROR WITH SLOPE FACTOR 


items, whereas the others are based 170 items. They may also 
due the less dependable data available for the 50-year floods. 

Comparison shown with the graphical study (using channel slope No. 
the same stages. Results the datatron and graphical method are not 
strictly comparable, first because the difference methods, second be- 
cause items were used the graphical analysis and 170 items the 
datatron analysis. However, the final datatron results, using the *85-10” 
slope, indicate practically the same standard error the channel slope 
No. used the graphical analysis. This gratifying because the “85-10” 
slope far simpler compute. 

Figure shows the typical regression equation with drainage areg and 
slope and plot the regression coefficients. The consistent variation 
regression coefficients evident. (The coefficient, for the drainage area, 
might well 1.00 throughout). The multiple correlation coefficients for 
each the regressions range between 0.93 and 0.97. 

Opportunity still exists for improvement the standard errors shewn 
Table the use other parameters employed the graphical study 
the improvement parameters such precipitation indices. Further 
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Table 2.- Multiple Correlation with Drainage 
Area and Slope (85 10%) 


Standard error 


Datatron Graphical Error 
log ave ave explained 
Ave. 


*Graphical standard error computed using channel slope No. 


work proceeding along these lines. This report presented progress 
report with definite finding one phase flood-frequency analysis, and 
not intended provide the final answer. 

The percent main-channel slope found effective for New 
England may not prove the best for other regions but should serve 
guide similar studies made elsewhere. 
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STORM WATER DRAINAGE THE CHICAGO AREA 


Horace Ramey,! ASCE 


SYNOPSIS 


This paper, read before the Hydraulic Division the American Society 
Civil Engineers February 25, 1958, discussion past and recent 
flooding conditions the Chicago area and some recommendations for their 
improvement. 

The opinions and recommendations stated herein are based fifty years 
investigation and experience with these problems. The recommendations 
discussed are not meant specific, but rather assist future students 
the subject arriving definite solutions. 


Floods have plagued mankind since the great Deluge 2349 B.C. That was 
about 4000 years ago. One might think that the Chicago area, much 
Continental Divide, with drainage both into Lake Michigan and into the 
Mississippi Valley, would comparatively free from floods, but such not 
the case. 

protect Lake Michigan, the source its public water supply, from con- 
tamination sewage, Chicago decided many years ago keep its drainage 
out the lake. This eliminated the quickest outlet for storm drainage, and, 
because the topography, placed the main outlet for such drainage the 
southwesterly corner the Chicago area. This, and the generally flat terrain 
the region, has produced many problems, some which are yet 
solved. 

the matter flood runoff, the Chicago region served three 
streams: the Chicago River, which drains the northern area the city itself 
far south 87th Street; the Des Plaines River, which drains the westerly 
areas; and the Calumet River, which drains the entire Calumet region south 
87th Street. 


Note: Discussion open until September 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1995 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. April, 1959. 
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knowledge the topography the Chicago region essential the 
study its flooding; and study the past furnishes the best data for the 
solution future problems. 


Topography, Drainage Areas 


The most important stream is, course, the Chicago River, which extends 
due west from Lake Michigan, little more than mile, the Forks. From 
this point the North Branch extends north northwest, about twenty-nine 
miles, its source the Skokie Valley, west Highland Park and Lake 
Forest, Lake County, The lower six-mile reach the North 
Branch the Chicago River has been improved for navigation and flood flow, 
Lawrence Avenue. Beyond that point the North Shore Channel The Metro- 
politan Sanitary District Greater Chicago extends eight miles con- 
nection with Lake Michigan, Wilmette. This channel (dry weather flow ca- 
pacity, 1000 cfs) passes the drainage the North Shore towns, Evanston 
Glencoe, into the North Branch. 

About three blocks north Lawrence Avenue dam was constructed 
hold the North Branch the Chicago River its natural state, hydraulically. 
Above this dam, the North Branch extends westerly and then northerly about 
eight miles and then divides into three small streams, the West Fork, the 
North Branch, and Skokie River. The Skokie Marsh covered five-mile reach 
this latter river, west Winnetka and Glencoe, and, its natural state, 
served retain certain flood water. Much this marsh has been drained 
and the area used for residential building; and lagoons have been constructed, 
control some extent the flow water the Skokie River. 

-The South Branch the Chicago River, extended from the Forks (Lake 
Street) due south about two miles, and then southwesterly two miles further, 
Ashland Avenue and 26th Street, where divided into the South Fork and 
the West Fork. The South Fork extended south one and one-half miles 39th 
Street and Racine Avenue. The West Fork extended west about two and one- 
half miles Kedzie Avenue and 32nd Street, the easterly limit slough 
known Mud Lake, which was the source the South Branch. Mud Lake ex- 
tended southwesterly, about six miles, Harlem Avenue and 49th Street, 
where connected with the Des Plaines River. 

The Main Drainage Canal was constructed from the West Fork the South 
Branch the Chicago River, Damen Avenue (one-half mile west Ashland 
Avenue) extend southwesterly across the Continental Divide just south 
Mud Lake, past Summit, and down the Des Plaines Valley, connect with the 
Des Plaines River Lockport, 

Present drainage areas the Chicago River (above Summit) are 309 square 
miles. 

The Des Plaines River, rising southeastern Wisconsin, west Racine, 
flows southward, parallel and about ten miles west the west shore 
Lake Michigan, about 49th Street and Harlem Avenue, Chicago (one mile 
south Riverside), and then makes sharp bend the southwest and con- 
tinues past Lockport and Joliet, Dresden, where unites with the Kankakee 
River form the River. Drainage area the Des Plaines River, 
above Riverside, 643 square miles; above Lemont, 690 square miles; and 
above Lockport, 713 square miles. 


J 
in 
n 


DRAINAGE 


The most important tributary the Des Plaines River Salt Creek, which 
comes from the west about 39th Street. Salt Creek rises Cook County, 
northeast Palatine, flows southward through the eastern section DuPage 
County, past Bensenville, Villa Park, and Elmhurst, Hinsdale; and then 
flows eastward past Western Springs and LaGrange Park, through Bookfield, 
into the Des Plaines River. Its drainage area about 152 square miles. 

Flag Creek rises the southern part Hinsdale and Western Springs and 
flows south about 87th Street, then southwesterly and into the Des Plaines 
River about half-way between Willow Springs and Sag. Its drainage area 
about square miles. 

The Little Calumet River rises Laporte, Indiana, and flows westward 
into about 170th Street, then northwesterly point Blue Island 
about one-half mile west Ashland Avenue, Then bends 
sharply the east and flows about seven miles junction with the Grand 
Calumet River, Torrence Avenue and 139th Street. From this junction the 
Calumet River (drainage area square miles) extends north about six miles, 
past Lake Calumet, Lake Michigan 90th Street. From this same junction, 
the Grand Calumet River (drainage area 50.1 sq. mi.) extends eastward the 
Lake Michigan shore near the east line Gary, Ind. Only sand bar sepa- 
rates this river from the lake. The Grand Calumet really bayou, having 
once been the outlet channel the Little Calumet. 

The Indiana Harbor Ship Canal connects the Grand Calumet River with Lake 
Michigan, through East Chicago, Indiana. The drainage area the Grand 
Calumet, east this canal, 30.8 square miles; and west it, 19.3 square 
miles. 

The original watershed the Little Calumet River, above Riverdale, was 
about 670 square miles. This included about square miles the eastern 
half the Sag Valley, drained Stony Creek, tributary the Little 
Calumet, and seven square miles the Blue Island area. also included 
about 370 square miles, Indiana, the upper reaches the Little Calumet, 
which, since 1925 have drained directly into Lake Michigan, near Dune Park, 
Indiana, through Burns Ditch. Deducting these areas, total 452 square miles, 
from the original 670 square miles, leaves the present watershed the Little 
Calumet, above Blue Island (its junction with the Calumet-Sag Channel) 
approximately 218 square miles. 

The drainage area the Little Calumet River, between Blue Island and its 
junction with the Grand Calumet, eight square miles, plus thirteen square 
miles tributary the Calumet storm water pumping station; total thirty-one 
square 

The present drainage area the Little Calumet River, reversed, 
flow into the Calumet-Sag Channel, Blue Island, 306 square miles. 

The drainage area the Calumet-Sag Channel, east the Continental 
Divide, about eighty-two square miles; and west this Divide, about twelve 
square miles. 

The most important tributary the Little Calumet River Thorn Creek 
(drainage area 115 square miles), which rises southwest Chicago Heights 
and flows into the Little Calumet the western part Lansing, One 
the principal tributaries Thorn Creek Butterfield Creek (drainage area 
about thirty square miles), which rises the southwest corner Cook 
County, west Matteson, and flows through Olympia Fields, and then south 
Flossmoor and Homewood, join Thorn Creek Glenwood. easterly 
tributary Thorn Creek Deer Creek (drainage area about sixteen square 
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miles), which rises Crete and flows northerly through East Chicago 
Heights, and then northwest along the Glenwood-Dyer Road, into Thorn Creek 
Glenwood. 

The area south 159th Street and west the Central Railroad 
drained the Calumet Union Drainage Ditch (drainage area about twenty-four 
square miles), which runs easterly from Kedzie Avenue, along the line 
163rd Street, and into the Little Calumet River South Holland. 


Topographic Characteristics 


The area discussed herein may generally described flat and rolling 
and although there are many ridges and divides the area, they are not obvi- 
ously distinct. 

The existence large flood plain evident those areas containing 
ponds and sloughs well the areas bordering many the shallow stream 
beds. 

Most the streams the area flow parallel the lake shore result 
the terminal morraines created the receding glaciers. 


Changes Stream Conditions 


The first change from the natural conditions these streams, regard 
storm flow, occurred 1848, when the Illinois and Michigan Canal was 
opened, from the South Branch the Chicago River Bridgeport (Ashland 
Avenue) through the Des Plaines Valley, Lockport. About year later, the 
Calumet Feeder the Canal went into service, from the Calumet 
River near Blue Island, through the Calumet-Sag Valley the Canal, 
Sag. 

Between 1867 and 1871, the and Michigan Canal was deepened, from 
Bridgeport Lockport, for gravity flow 1,000 cfs from the Chicago 
River, lieu about 250 cfs previously pumped. 

1871, the Ogden-Wentworth Ditch was dug, from the Des Plaines River, 
through Mud Lake, the West Fork the South Branch the Chicago River. 
The purpose was drain Mud Lake, enhance local real estate values. This 
ditch was feet wide and its bottom was below the bed the Des Plaines 
River. the spring 1872, the Des Plaines River floods greatly enlarged 
the new ditch, depositing the soil the Canal. 

Thereafter, more flood water than normal came from the Des Plaines 
River into the West Fork the South Branch. Even the dry weather flow 
the Des Plaines River followed the same course and then flowed down the 
Canal. This condition continued until June 1877, when the Ogden Dam was 
built. 

The Ogden Dam was row sheet piling, feet long, with crest ele- 
vation +11.8 CCD, which was about two feet lower than the natural banks 
the Des Plaines River, and foot two higher than the surface Mud Lake. 
was built the arm the Des Plaines River which led Mud Lake, west 
Harlem Avenue, near 49th Street. This dam was rebuilt 1885, stone- 
filled timber cribs, with its crest the same elevation. The new dam served 
until 1894, when the Riverside Spillway, 397 feet long, elevation +16.25 


CCD, was built the east side the Des Plaines, short distance north 
Ogden Dam. 
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AND VICINITY 


The Riverside Spillway was closed 1908 when the Willow Springs Spill- 
way was built, replace it, and discharge like amount flood water from 
the Des Plaines River into the Drainage Canal, instead into the South 
Branch the Chicago River. The Willow Springs Spillway was closed 
1955, and since that date Des Plaines River floods have not affected the flow 
either the Chicago River the Drainage Canal. 


Chicago Drainage Canal 


The greatest change stream conditions the Chicago region was the 
construction the Chicago Drainage Canal, 1892-1900, from the West Fork 
the South Branch the Chicago River, Damen Avenue, twenty-nine 
miles, the Des Plaines River Lockport, This was later extended, 
1907, four miles, connect with the Upper Basin the and Michigan 
Canal and the Des Plaines River, the northern part Joliet. This channel 
was designed, 1892, with normal hydraulic capacity flow 10,000 cfs 


water, carry the runoff the maximum storm expected the Chicago area 
for the next thirty years. Its capacity, the rock section, fifteen miles, from 


Willow Springs Lockport, 14,000 cfs. 

Construction the Drainage Canal was accomplished the improvement 
the South Branch the Chicago River, capacity flow 10,000 cfs 
water velocity less than feet per second. 


North Shore Channel 


The North Shore Channel was constructed, 1908-1910, from Lake Michigan, 
Wilmette, southwest about two miles Emerson Street, Evanston, and 
then due south the line Kedzie Avenue, six miles, Lawrence Avenue. 
joined the North Branch the Chicago River about three blocks north 
Lawrence Avenue. Its nominal flow capacity was 1000 cfs. 

Prior the construction the North Shore Channel, the marshy land 
around Emerson Street, and southward, drained south through small brook 
into the Chicago River, south Foster Avenue. The North Shore Channel 
added significant drainage area the North Branch the Chicago River, 
but did speed the runoff flood water. 


Calumet-Sag Channel 


The Calumet-Sag Channel was constructed, 1911-1922, from the Little 
Calumet River near Blue Island, sixteen miles through the Sag Valley, con- 
nect with the Main Drainage Canal, Sag, about three miles upstream from 
Lemont. The nominal flow capacity this channel was 2000 cfs. Its depth 
was feet and its minimum width feet. This channel now being 
widened, the Federal Government, for navigation, bottom width 225 


feet, depth feet below water. This will increase its nominal hydraulic 
capacity 5000 cfs. 


Summary Regarding Man Made Changes—General Effects 


The areas discussed this paper are densely populated and highly industri- 
alized. Some the stream characteristics have been altered dredging and 
the addition many man made structures. 
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effect these man made changes has been change, some cases, 
the outlying borders the drainage basins. This, together with the general 
topographic characteristics gently sloping plains and flattish slopes, pre- 
sents drainage problem which requires bold and imaginative solution. 

Considerable progress has been made the past with the creation the 
existing open channel system, well the sewer program presently under 
way. is, however, obvious that the tremendous population explosion and 
rapid rate urbanization require immediate action remedy deteriorating 
situation. 


Introduction Flood History 


Considerable material describing early floods the Chicago region 
available the form historic rhetoric. 

Some rainfall data far back 1875 available for use. This data is, 
however, not sufficiently complete afford any scientific analysis their ef- 
fect. The following briefly summarizes the general impressions those re- 
corded floods and periods intense rainfall. 


Early Floods 


The first record flood the Chicago region was March 29, 1674, 
the explorer priest, Marquette. and two other Frenchmen and number 
Indians, their way from Green Bay the principal town the 
Indians, had spent the winter, after December 1673, near Damen Avenue 
and the West Fork the South Branch the Chicago River. March 29, 
1674, was driven from his camp high water coming through Mud Lake, 
from the Des Plaines River, with the spring break-up the ice. proceed- 
ed, the next day, the South Branch, through Mud Lake, without any portage, 
and then down the Des Plaines River, which was flood (Brown’s 
“Drainage Channel and Waterway,” pages 103-104). 


Flood March 12, 1849 


most disastrous flood occurred the Chicago River, March 12, 1849. 

precipitation records are available, but this flood followed the pattern 
previous floods coming from the Des Plaines River, through Mud Lake, into 
the South Branch. ice jam formed the South Branch, south (upstream) 

the Madison Street bridge and raised the water and floating ice two three 
feet. The jam broke A.M., March 12, and the rapid current swept away 
every vessel moored below that point. The mass ice and sailing vessels 
struck the Madison Street bridge, wrecked it, and added the bridge the float- 
ing debris. Next, the Randolph Street bridge was swept away, were 
number canal boats and sailing vessels moored the South Branch, north 
Randolph Street. 

The mass ruin then swept into the main river, where the banks were 
lined with steamers, schooners, brigs, canal boats, scows, etc. The Wells 
Street bridge was open and the middle the stream clear, but the floating 
mass vessels and wreckage extended from bank bank and soon swept 
away both ends this bridge and mingled with the general wreckage. The 
larger vessels crashed together, entangling booms, bowsprits, and rigging. 
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The smaller vessels were The ice, places, was wedged several 


feet above water level. 


The Clark Street bridge was completely swept away and great many 
the larger vessels became wedged together short distance east Clark 
Street, and with the general wreckage formed dam, o’clock noon. 
Above this dam, ice and wreckage accumulated considerable height; and 
below it, numbers vessels were swept into Lake Michigan, many which 


were wrecked. 


The dam still held 10:00 A.M. the following day, but most the ice had 
passed under the vessels and the water level had fallen about two feet. The 
crisis was over. Fears that the ice the North Branch would break and 


add the damage, were not realized. 


The flood was general and caused considerable damage the Des Plaines 


and Valleys. 


Other Floods 


Other floods occurred Chicago, due the following extremely heavy 


rains: 


September 1875 
July 25-26, 1878 
July 6-7, 1879 


November 11-12, 1881 


November 5-6, 1883 
March 25-26, 1884 
June 2-3, 1885 


However, the records available are not sufficient for any study these 


floods. 


Summary Runoff from Past Floods 


For convenience study, the principal data regarding past floods the 


Chicago area are listed the tabulation following: 


August 2-3, 1885 
May 6-8, 1892 


August 11, 1908 
April 29, 1909 


August 14-15, 1909 
March 17-18, 1919 


May 3-6, 1919 
Sept. 11-13, 1936 


July 1938 
July 1943 


Rainfall 
Inches 


4.14 


Time 
Hours 


Www ww WwW 


Rainfall 


inches 
inches 
inches 
inches 
inches 
inches 
inches 


Runoff 
cfs 


9,000 
8,900 


10, 600 


15,000 


16, 000 


6.33 
2.92 
3.52 
2.77 
4.58 


April 5-7, 1947 5.27 29,100 
March 19, 1948 3.50 24, 500 
July 14-15, 1949 4.64 
July 16-17, 1950 5.40 
May 10-11, 1951 3.91 16, 000 
October 1954 6.66 
7.20 
July 12-13, 1957 6.00 28, 000 
6.34 


The conclusion drawn from study these floods that the 
Chicago area within the past thirty years the runoff, from one cause an- 
other, has increased least two and one-half times. The Chicago Drainage 
Canal, planned 1892 carry the runoff from the greatest flood which might 
occur the Chicago area, then estimated 10,000 cfs, has now become in- 
adequate. Many changes have taken place the Chicago area since that time, 
sixty-five years ago, and the whole storm drainage problem must re- 
evaluated. Apparently outlet capacity 30,000 35,000 cfs will needed. 


Floods the Des Plaines River 


The Des Plaines River has drainage area 643 square miles, above 
Riverside. this, some 152 square miles are the watershed its princi- 
pal tributary, Salt Creek. The drainage area Salt Creek, above the Fullers- 
burg dam, just north Hinsdale, about 122 square miles. 

From the Hoffman dam, Riverside, the Des Plaines River extends thirty- 
one miles north the Cook-Lake County line, just north Wheeling, and 
fifty-seven miles further its source Wisconsin, west Kenosha and 
Racine. Downstream from Riverside this river extends southwesterly, about 
twenty-six miles, Lockport. The slope the flood plain the river, and 
the river itself, averages about 1.0 feet per mile between Riverside and the 
north line Cook County, and about 1.9 feet per mile between Lockport and 
Riverside. 

The average cross section the Des Plaines River, between Riverside and 
North Avenue, for dry weather flow, slightly less than 400 square feet; and 
the section between North Avenue and Wheeling probably less. This means 
channel about 150 feet wide and less than feet deep. 

The Forest Preserve District Cook County, conjunction with the State 
has adopted project control the low water flow the Des 
Plaines River, which times low cfs, for weeks end, the 
thirty-one mile stretch between Riverside and the Cook-Lake County line. 
This involved the construction twelve thirteen low head dams, provide 
least two feet water above each dam. Dams Dempster Street, Devon 
Avenue, and Touhy Avenue were built 1946. The Hoffman dam, Riverside, 
was rebuilt 1950-1956. Dams Hintz Road, one mile south Wheeling, 


and Foundry Road, were replace existing dams. One other dam was built 
1957. 
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least twenty-nine highway and nine railroad bridges cross this stretch 
the Des Plaines River, the average horizontal clearance such bridges 
being about 150 feet. The average loss head through these bridge openings 
times high water, may inches per bridge. The flood stage 
North Avenue may therefore feet higher because losses through 
these restricted openings. The flood stage Wheeling may increased 
twice much. Accumulated drift bridges and other narrow places has, 
doubt, caused further losses, for there record this river having been 
cleared during the past twenty-five years. 

Some recent floods the Des River, measured Riverside, 
are follows: 


Maximum Runoff 
Riverside Discharge cfs 
Gauge cfs per sq. mi. 

April 21, 1881 13, 500 21.0 
February 1887 20.2 10,324 
June 24, 1892 20. 
March 25, 1904 21.45 10, 780 16.8 
March 1908 21.47 10, 810 16.8 
May 1909 21.97 11,620 18.1 
March 18, 1919 7,450 11.6 
May 1933 5,200 8.1 
Sept. 17, 1936 5,160 8.1 
July 1938 5,350 8.3 
March 16, 1944 4,600 
April 1947 780 
March 20, 1948 6,510 10.1 
April 26, 1950 340 9.9 


The measurement, April 21, 1881, was made Lyman Cooley, flow 
over weir. The measurement, February 1887, was made 
Johnson, using rod floats. 

The flood April 1881 was due the melting ice and snow, amounting 
about one foot water, which had accumulated ground which was practi- 
cally saturated when winter set in. 

Four floods Salt Creek, measured the Fullersburg dam, 1876, 
1881, 1886, and 1887 produced flows about 2800 cfs. these cases, the 
runoff was about 22.6 cfs per square mile, from the 122 square miles the 
watershed above Fullersburg. The actual measurement 1881 was 2760 cfs. 

The floods 1881, 1892, and perhaps 1909 may classed great 
floods. fortunate for the residents the Des Plaines Valley that 
great flood has occurred the Des Plaines River since 1909. The highest 
water levels this time were March 1919; and these levels were about 
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3.5 feet below the levels due the 1881 flood. The recorded discharge 
Riverside, March 18, 1919, was 7450 cfs. However, with discharge 6510 
cfs March 20, 1948 and 6340 cfs April 26, 1950, the water levels the Des 
Plaines River between Madison Street and North Avenue closely approached 
the high levels March 1919. This indicates some obstruction flow and 
the need for clearing the channel. 


report the “Improvements the Des Plaines River” was made July 


1924, Ernest Cooley, Hydraulic Engineer, Edward Kelly, Chief 
Engineer The Sanitary District Chicago. For the critical reach between 
the Hoffman dam, Riverside, and North Avenue, the improvement recom- 
mended was follows: 


Lower the crest the Hoffman dam from +25.56 elevation +20.00 

Excavate the channel from the Hoffman dam, 4800 feet the junction 
with Salt Creek, bottom width 300 feet and bottom grade 
+15.00; 

Excavate the channel, from the Salt Creek junction point about 
mile north North Avenue (Elmwood Park), bottom width 200 
feet and bottom grade rising +17.00 22nd Street, +20.00 Lake 
Street, and +22.00 North Avenue; 

Construct sub-dam 22nd Street with crest elevation +22.00; 

Lower the crest the Public Service dam, below Lake Street, from 
+28.54 elevation +25.00. 


This plan, estimated cost several millions dollars, would lower the 
flood water stage feet, from the Riverside dam North Avenue and 


The dams would cause low water boating pools, 18,000 feet long 


through Riverside 22nd Street, 16,000 feet long from 22nd Street Lake 
Street, and 15,000 feet long from Lake Street Elmwood Park. 

This was good plan 1924 and good plan today, with proper modi- 
fications, including its extension the north line Cook County. Recent 
studies indicate that the bottom grade the channel should about +13.00 
C.C.D. from the Hoffman dam the mouth Salt Creek and then rise, ona 
slope 0.002, about +17.20 12th Street, +21.20 North Avenue, +33.00 
Golf Road, +36.00 Lake-Euclid Avenue, and +45.00 the north line 
Cook County. The depth, below high water, would feet. Bottom width 
the new channel would 300 feet, from the Hoffman dam the mouth Salt 
Creek, then about 230 feet North Avenue, then 210 feet Higgins Road, 
diminishing 150 feet Palatine Road. Changes channel widths would 
made points where tributary creeks enter the Des Plaines River. main- 
tain desired low water levels, sub-dams should constructed, proper 
intervals. lessen the rise water level over these dams, pass high 
flood flows, the dams could built diagonally instead directly across the 
stream, and thus given greater crest lengths. Growth shrubbery, which 
would obstructive flood flow, should restricted. 

believed that channel improvements, such those suggested above, 
would enable the Des Plaines River carry any flood which may reasonably 
expected the future. The depth cut, below the bed the present 
channel, would probably average feet, between the Hoffman dam and 
North Avenue, and perhaps feet between North Avenue and the County Line. 
The construction cost might $6,000,000 $9,000,000. 
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Improvement would also needed Salt Creek, Addison Creek, Silver 
Creek, and other tributaries. Reports have been made, the Division 
Waterways the State indicating that such construction might 


Addison and Silver Creeks $4, 790, 520 
Flag Creek 836,170 


Flag Creek (drainage area about square miles), which rises the 
southern part Hinsdale and Western Springs, north 47th Street, 
small stream, about seven miles long. The Dlinois Toll Road (now under con- 
struction) parallels Flag Creek, from Street, north about 47th Street, 
and crosses it. Considerable flooding has occurred, the past, Hinsdale 
and Western Springs, because the inadequate natural channel Flag Creek. 
recent years considerable amount building has been permitted 
low lands subject overflow from high floods the Des Plaines River. Inno- 
cent purchasers houses such areas have suffered unexpected and un- 
warranted damage from high water. public agency should approve such 
construction, regardless the pressure applied subdividers and builders. 


Little Calumet River 


study floods the Little Calumet River indicates that storm runoff 

rates the Calumet area have least doubled the past fifty years. The 

drainage area this stream has been subject several changes, over the 
ind Years. The original 670 square miles the Little Calumet above Riverdale, 
included seventy-five square miles the Sag Valley and seven square miles 


around Blue Island, now tributary the Calumet-Sag Channel, well 370 
square miles Indiana which now drain Lake Michigan through Burns 
Records the worst floods are follows: 

Drainage 

Area Rate per 

Runoff Sq. Mi. Sq. Mi. 

Riverdale 

Salt March 1908 11,000 cfs 670 16.4 

June 11, 1939 3,440 cfs 198 17.4 cfs 
March 16, 1944 198 21.7 cfs 

Blue Island 

April 5-7, 1947 9,100 29.8 cfs 
March 19, 1948 5,000 cfs 306 16.3 cfs 

October 9-11, 1954 6,700 cfs 306 21.9 cfs 

bly 


new lock and control gates are built the Calumet River part 
the project for widening the Calumet-Sag Channel, now under construction. 

This will leave only the thirty-eight square miles the Calumet River water- 
shed tributary Lake Michigan. The remaining 268 miles the present 306 


July 12-13, 1957 7,200 cfs 23.5 cfs 
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square miles watershed the Little Calumet, Grand Calumet, and Calumet, 
will tributary the Calumet-Sag Canal. 

For present flood relief into Lake Michigan, channel 3,500 4,000 
square feet cross-section available. the future, only the capacity 
the control gates the new lock will available. 


would appear that the minimum capacity the Calumet-Sag Channel, for 


the passage the runoff from great flood, should about 10,000 cfs. This 
much, more, could come from runoff cfs per square mile, from the 
306 square mile watershed east Blue Island and the square miles the 
Calumet-Sag valley. About per cent the 221 square miles within the 
city limits Chicago, square miles, south 87th Street, are included 
the area draining the Little Calumet River and the Calumet-Sag Channel. 
The channel the Little Calumet River, upstream from Blue Island, 
quite inadequate for storm water flow. 1941 the State improved 
this channel for 50-foot width, from the Indiana State line the mouth 
Thorn Creek; and 1946 extended this 50-foot improvement Roll Street, 
Blue Island. Under present conditions greatly increased runoff due 
the considerable building the area, appears that the width this channel 
should more than doubled and perhaps the depth should increased. Open- 
ings the many bridges should greatly increased. pending improvement 
Thorn Creek, already authorized the State will increase the 
rate storm water flow; and direct attention the need for further im- 
provement the Little Calumet River above Blue Island. 


Summary—Floods: Need for Remedial Measures Obvious 


The foregoing has outlined chronologically, the history storms and re- 


‘sultant floods the Chicago region. can seen, there exists some data 


recent years which can used appraise and investigate few past 

storms and their relation flooding. becomes more obvious upon such ex- 
amination that immediate remedies must found for the flooding situation 
the area. The flood frequency becoming greater with successively smaller 


rainfall and promises become still greater with further urbanization and 
growth. 


City Chicago 


The City Chicago embarked program construction storm relief 
sewers, 1947, expected entail ultimate expenditure $165,000,000, 
catch with existing needs. Three bond issues provide funds for this 
construction have already been approved the voters, follows: 
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hoped that this work can carried the rate about $15,000,000 
per year, 1957 and after. 

addition the work the City Chicago, The Sanitary District 
Chicago expended $47,367,000 its system South Side relief sewers (to 
December 31, 1956); and currently spending about $5,200,000 ex- 
tension the Racine Avenue Pumping Station, for pumping storm water only. 
The total expenditure, the Sanitary District, for this storm drainage will 
about $52,570,000. 

This system sewers with the existing 39th Street Conduit was designed 
carry 5500 cfs storm runoff, from 27.6 square miles, between 31st 
Street and 87th Street, and east Racine Avenue Lake Michigan; 
runoff rate 200 cfs per square mile. These sewers have capacity carry 
400 cfs sewage, addition the 5500 cfs storm flow, total 5900 cfs. 

The sewer outlet capacity the City Chicago, past and prospective, in- 
cluding the South Side sewer system the Sanitary District, indicated 
the tabulation following. Runoff rates are based the present 221 square 
miles area within the City Chicago. 


Runoff 

Years Sewer Outlet Capacity cfs per sq. mi. 
1900-1910 10,000 cfs 45.5 

1920 14,000 cfs 64.0 

1930 20,000 cfs 91.6 

1940 30,000 cfs 137.0 

1950 35,000 cfs 159.0 

1957 37,000 cfs 173.0 

1960 41,000 cfs 187.0 

1970 55,000 cfs 250.0 


Sewer outlet capacity, thirty years ago, was about ten times the dry 
weather flow. The outlet sewers now being constructed are designed with ca- 
pacities thirty fifty times the dry weather flow. This will reduced 
high water levels the Chicago River and the Drainage Canals; but the 
storage capacities these large relief sewers will considerable. 

These present and prospective sewers are designed flow their full ca- 
pacities into the rivers canals water level +1.50 CCD. This water 
level may reach +4.00 the Chicago and Calumet Rivers and the effect 
this would reduce the flow capacities these sewers per 
cent their designed capacities; reduce the total (1970) capacity from 
55,000 cfs about 40,000 cfs. 

not expected that the flow capacities all these outlet sewers could 
ever developed the same time. The times concentrations floods 
the various sewers and sewer systems range from less than thirty minutes, 
three hours. However, 6-inch rain six hours, such fell the 
entire Chicago area July 12, 1957, could, through this completed sewer 
system, deliver amazing amount storm runoff into the Chicago and 
Calumet Rivers and the connecting canals. July 1957 lot this went into 
the basement householders, because the inadequacy the local sewers. 

doubt these sewers, particular, the recent auxiliary outlet sewers, 
have been designed liberal basis. However, they have not been greatly 
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over-designed, study how much flood water they can deliver the outlet 
streams may more realistic method estimating the volume future 
great flood, than the method outlined previous pages, estimating from 
the amounts water which have been discharged through the Drainage Canal, 
times past floods. 

study these outlet sewers, only the 221 square miles within the city 
limits Chicago, indicates the following: 


Percentage Sewer Capacity cfs 
Sewer Capacity Total 75% 50% 

North Lawrence Avenue 10% 5,500 4,120 

Lawrence Avenue Lake Street 17% 9,350 7,010 4,670 

Total North Branch 27% 14, 850 11,130 

South Branch Chicago River 30% 16, 500 8,250 

Total Chicago River 57% 23,500 15,670 
Drainage Canal, Western Ave. 

Harlem Ave. 16% 800 600 

Total Main Channel 73% 40,150 30,100 20,070 

Calumet Area, South 87th Street 27% 14, 850 11,130 7,430 

Total, City Chicago 100% 55,000 41,230 27,500 


Storm Flow—Computed from Sewer Capacity 


may assumed, quite probable, that least per cent the total 
flow capacity these sewers will developed simultaneously, some future 
flood, and that the runoff from the 221 square miles within the present 
Chicago city limits will least 27,500 cfs. this must added the 
storm runoff from about 100 square miles, north Chicago, the upper 
reaches the North Branch and the North Shore Channe, including 
Evanston and the other North Shore towns; and the runoff from about 260 
square miles tributary the Little Calumet River, south the city. Such 
runoff may reasonably calculated cfs per square mile; and would 
amount about 3000 cfs the North Side, and 7800 cfs the Calumet 
region. 

The 3000 cfs the North Side, added the 20,070 cfs total flow the 
Main Channel, from the table above, would indicate storm flow least 
23,070 cfs the Main Drainage Canal, above Summit. This runoff, from the 
308.9 square miles drainage area above Summit, would average 
rate about cfs per square mile. 

The 7800 cfs the Calumet region, added the 7430 cfs flow the 
Calumet-Sage Channel, from the table above, would indicate storm flow 


about 15,230 cfs the Calumet-Sag Channel. This runoff, from the 306 square 


miles tributary the Calumet-Sag Channel, would average rate 
about cfs per square mile. Storm flow from the ninety-six square miles 
the Sag Valley ignored this calculation, since much this could have 
already passed downstream before the peak concentration the entire 
Calumet region. 

The total storm flow, under this method calculation, would 23,070 cfs 
from the area north 87th Street plus 15,230 cfs from the Calumet region, 
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38,300 cfs. This can compared with the conclusion 30,000 35,000 
cfs derived (page 25) from study the runoff from past floods. 


Introduction Recommendations 


The following discussions and tabulation comprise program action 
geared increase the outlet capacities for the expected flood waters the 
years come. They are means the total solution the problem but 
only make available sufficient outlet capacities for the rapidly expanding 
sewer programs. should understood that these recommendations are 
based the data and experience taken from the operation the open channel 
system and are not meant solve localized problems due inadequate sewer 
capacities. This program suggested that part the total program which 
will fulfill the outlet capacity needs for flood waters for the next fifty years. 


Main Outlet Channels 


Since storm runoff the Chicago and Calumet areas has approached 30,000 
cfs three different occasions the past ten years (March, 1948, October 
1954, July 1957), and since building and paving are increasing very sub- 
stantial rates, apparent that the runoff from great storm may 30,000 
35,000 cfs, the not distant future. obvious that outlet capacity for 
such storm flow should least 30,000 cfs. 

Chicago the one city, the Great Lakes, the fortunate position not 
having drink its own sewage, pass for others drink. However 
dilute might be, the thought repugnant. For the past thirty-five years, 
since August 1922, when the Calumet-Sag was opened, Chicago has kept all 
its sewage out Lake Michigan, the source municipal water supply. has 
not been possible keep the peaks flood flows from the Calumet Rivers 
out the lake, whenever such flood flows have exceeded the hydraulic capaci- 
the Calumet-Sag Channel. The same true the North Shore Channel, 
where storm runoff recent years has had great increase, and where the 
consequent water level has, many occasions, been high enough cause 
flow over the restraining walls Wilmette. only two occasions (October 
1954 and July 1957) has been necessary permit the Chicago River flow 
into Lake Michigan, prevent damage Loop property, from high water. 
These were occasions what might called great storms. 

may that this policy Chicago, cherished over the years, keeping 
its drainage out its drinking water, may have abandoned; least dur- 
ing great storms. recommended, however, that sufficient outlet capacity 
may provided for storm flows that relief discharge into Lake Michigan 
can limited the peak runoff from great storms. 

that end, hydraulic study has been made the present outlet channels 
and how much their flow capacities should increased, meet the needs 
the near future, discharge capacity 30,000 cfs, assuming that any flood 
flow excess that amount could spilled into Lake Michigan. This study, 
which should considered preliminary only, indicates need for the 
following improvements: 


Main Outlet, Drainage Canal, from Sag Lockport— 
Widen present channel, from 160 feet 300 feet, more. 
This channel would carry about 30,000 cfs, with water surface -3.0 
Sag and -8.0 Lockport. 
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Calumet Outlet, Calumet-Sag Channel, from Blue Island Sag— 
Present widening 225 feet, Government. This channel should 
carry about 10,000 cfs, with water surface -3.0 Sag and +4.0 
Little Calumet River. 


Chicago Outlet 
Drainage Canal, from Willow Springs Sag— 


Widen present channel, from 160 feet 230 feet. This channel 
should carry about 20,000 cfs, with water surface -1.5 Willow 
Springs and -3.0 Sag. 


Drainage Canal, from Summit Willow Springs— 
improvement indicated. Present channel should carry about 
20,000 cfs, with water surface +0.1 Summit and -1.5 Willow 
Springs. 


Drainage Canal, From Western Summit— 
Widen present channel feet. The channel should then carry 
about 20,000 cfs, with water surface +1.5 Western Avenue and 
+0.1 Summit. Perhaps this reach should improved the 
standard dock section, with 240 ft. width between docks; and 
proper deepenings. 


Chicago River and Drainage Canal from Lake Street Western Avenue— 
improvement indicated. Present channel should carry about 
20,000 cfs, with water surface +4.0 Lake Street and +1.5 
Western Avenue. 


These estimates are based maximum stage +4.0 CCD the Chicago 
and Little Calumet Rivers, -3.0 Sag, and minimum -8.00 Lock- 
port. Holding the water level Sag about -3.00 very important. the 
March 1948 flood, this level was about +1.60. the October 1954 flood, 
was probably much higher, but the gauge failed. the July 1957 flood, was 
+0.40. These high water levels Sag cause unreasonable decrease the 
flow capacity the Drainage Canal above Sag; and the Chicago River. 

Further improvement the Chicago River and the Drainage Canal be- 
tween Damen Avenue and Summit could easily accomplished dredging 
the central 100 feet the cross-section greater depth. 

The improvement these channels, indicated, would involve the exca- 
vation approximately 7,200,000 cubic yards rock and about 4,300,000 
cubic yards glacial drift, earth, cost about $30,000,000. more 
detailed study would result better hydraulic balance between the various 
channels and sections channels. 


Efforts 1955 1957 


Since 1955, when the storm water drainage program within The Sanitary 
District Chicago, outside the City Chicago, then estimated cost 
$72,000,000, was submitted the General Assembly, little has 
happened change the picture. 

The Metropolitan Sanitary District Greater Chicago has caused the 
dredging 117,325 cubic yards material from the North Branch the 
Chicago River and the North Shore Channel, cost $378,960. This 
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way enlarged this channel, but such maintenance dredging did restore the 
original cross-section and improved conditions for flood flow. 

Currently, maintenance dredging the North Shore Channel, from Foster 
Avenue Wilmette, under way. Some 250,000 cubic yards silt and fill 
are removed, during the year 1958, cost about $1,110,000, 
restore the balance the North Shore Channel its original cross-section. 

Debris and silt and some shrubbery were removed from the channel the 
Des Plaines River between the old Lockport Controlling Works and Ninth 
Street bridge, Lockport, and the original channel restored 1956. Similar 
work now under way between the Ninth Street bridge and the Elgin, Joliet, 
and Eastern Railway bridge. This will finished 1958. 

Three additional sluice gates the old Lockport Controlling Works were 
equipped with operating machinery, 1956-1957, thus adding about 7,500 cfs 
the discharge capacity from the Main Drainage Canal, Lockport. Current- 
ly, studies are being made methods obtain 10,000 12,000 cfs more out- 
flow capacity, probably the Lockport Power House. Outflow capacity was 
somewhat inadequate during the October 1954 flood, when two the sluice 
gates were blocked some loose barges. 

The channel Deep Run being enlarged between Sixteenth Street bridge, 
Lockport, and its junction with the Main Drainage Canal below the Lockport 
Power House. The improved channel will have storm flow capacity 6500 
cfs. This work will finished 1958, cost about $500,000, and 
should alleviate the worst flood conditions Lockport, 

The City Chicago continued construction its long range program 
Auxiliary Outlet Sewers, presumably rate about $15,000,000 per year. 


bond issue $8,000,000, for this work, was approved the voters, 


1957, leaving about $70,000,000 the future financing yet approved. 


Estimates Cost Future Work 


little has been done toward the alleviation flood conditions Cook 
County, outside the City Chicago, that the estimates costs, submitted 
1955, may well repeated, with minor modification. One important and cost- 
item, however, must added, namely, the enlargement the Main Outlet 
Channel, $30,000,000. 


The entire future storm drainage program, not yet financed, outlined 
the tabulation following. 


Storm Water Drainage Program February 1958: 
Enlargement Main $30,000, 000 


Improvement Des Plaines River, from 


Riverside north line Cook County......... 9,000, 000 


Improvement North Branch Chicago River, 


Enlargement North Shore 2,000,000 


Enlargement Stony Creek, and construction 
tributary ditches and 8,000,000 
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Improvement Calumet Union Drainage 2,000, 000 
Improvement Feehanville and Weller 1,000, 000 
Improvements along Addison and Silver 4,000, 000 
Storm water outlets for Posen, Midlothian, 

Storm water outlets, etc. for villages 

Grand Calumet and Little Calumet 6,000, 000 

Total Cook County, outside Chicago....... $72,000, 000 

Total City Chicago, 70,000, 000 

Total storm water $142,000,000 


Storm water sewer systems for the various villages Cook County, esti- 
mated 1955 cost $31,000,000, are not included the above estimate. The 
improvements Salt Creek, DuPage County, estimated 1955 $5,000,000, 
some importance the Chicago area, but omitted. 


CONCLUSIONS 


From study floods the Chicago area, obvious that the present 
outlet channels are not adequate handle the storm runoff during great 
flood. change has been made the Main Outlet, that is, the Drainage 
Canal between Sag and Lockport, since 1900. that time, the storm runoff 
from the Chicago and Calumet areas has probably more than tripled. Local 
drainage has been and being improved many parts the area, but these 


improvements can not fully effective until the Main Outlet made adequate. 


Therefore, the enlargement the Drainage Canal might considered 
first importance. This will particularly true the stretch this canal, 
downstream from Sag, since the widening the Calumet-Sag Channel will, 
the next four years, cause tremendous increase storm flow from the 
Calumet region. Enlargement the Drainage Canal estimated cost 
$30,000,000, about the amount damage estimated from one great flood, 
October 1954. Damage from the flood July 1957 was less, but figure 
$10,000,000 has been mentioned. probable that the cost flood pro- 
tection would the cost flood damage. 

All the projects listed the program above, except the enlargement the 
Drainage Canal, are urgently needed for protection against the ordinary flood, 
that due ten-year storm. Residents some the communities in- 
volved were pitiful distress during the great floods 1954 and 1957. 

not the purpose this paper suggest specific improvements for 
handling storm water flow any part the Chicago area; but rather out- 
line the problem and furnish background material and factual data which will 
assist future students the subject arriving definite recommendations. 
Proper solution the problem will require construction program far beyond 
the imagination everybody except few who have been close enough 
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realize its magnitude. Considerable time, least ten years, will required 
for the physical construction, orderly fashion, the needed works, 
under any method financing which will not too burdensome. 

program which must necessarily carried out area-wide basis. 

Chicago great city, although young city, only 120 years old. 

Perhaps has grown too fast and some its problems have been neglected. 

can become the greatest city the world; and the solution the storm 

drainage problem, now, would enhance that prospect. the powers that 

will face that essential issue, and act, future generations will bless them. 
they procrastinate and postpone such really necessary improvements, future 

generations will condemn them for lack courage. well known that inde- 
cision the greatest failing the average office holder; but there time 
for indecision matter important large and populous community 

the prevention flooding. The longer such matters are deferred, grow- 
ing community, the greater the cost will be. 
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HYDRAULIC ANALYSIS SURGE TANKS DIGITAL COMPUTER 


SYNOPSIS 


shown how certain surge problems have been successfully adapted for 
solution IBM 650 computer. Machine results for the Oahe surge system 
are presented and compared with those obtained independent manual 
methods arithmetic integration. Machine-computation time for each the 
cases discussed averaged between one one and one-half hours. Sufficient 
data are given permit individual checks machine analysis and compu- 
tations relatively short time. 

Computer results are discussed indicate that treatment the total 
problem transients one separable parts the analysis surge-tank 
stability conventional methods leads inconclusive results. appears 
that comprehensive study digital computer the total problem 


transients feasible, but machine programming job great scope will 
required before fruitful results are available. 


INTRODUCTION 


High-speed, automatic digital computers have the potentiality opening 
new horizons understanding problem phenomena. Machine solutions will 
often provide broader perspective problem phenomena than can feasi- 
bly obtained manual computations because the relative rapidity with 
which the effects various factors and assumptions involved problem 
can tested and explored. demonstrating how this advantage digital 
computer has been actually utilized for solving surge problems particular 


case, additional evidence provided the growing usefulness computers 
engineering work. 


Note: Discussion open until September 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1996 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. April, 1959. 


Hydrs. Section, Army Engr. Div., Missouri River, Omaha, 
Nebr. 
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The study discussed herein concerned with mass-oscillation phenomena 
surge system, which includes pressure conduit and two restricted- 
orifice type surge tanks tandem independently connected the main con- 
duit (Fig. 1). The system supplies water reaction-type turbine. Water- 
hammer effects are not included the study. Thus, the study primarily 
concerned with the solution surge-tank problems, involving load rejection, 
load demand, and hydraulic stability the tanks under small load changes 
minimum head. 

Given the initial-flow conditions and all other pertinent data, the values 
various discharges and heads function time for the unsteady flow con- 
ditions within the surge system, resulting from power-load changes, may 
determined conventional methods arithmetic integration. (1,2)* Appli- 
cation these methods surge problems based upon the fundamental as- 
sumption that the compressibility water and the elasticity flow bounda- 
ries may ignored. This means that pressure changes the surge system 
are transmitted with infinitely high velocities and that water columns differ- 
ent parts the system may regarded solids evaluating positive 
negative accelerative effects. The same fundamental assumption used 
basis for development machine methods. However, since digital computer 
employed obtain solutions, there need make many other 
simplifying assumptions are ordinarily dictated the time limitations 
manual solutions. For example, acceleration heads within the risers, vari- 
ation friction factors with changing Reynolds Number, tee losses, and finer 
time increments the step-by-step process may all introduced with only 
comparatively small increase the machine time required for problem so- 
lutions. 

list the physical relationships required for machine solution surge 
problems associated with the system depicted Fig. and the sign convention 
that has been used are given Appendix III. 


Basis Analysis and Computation 


Head and Discharge Relationships 


Based upon the fundamental assumption instantaneous transmission 
pressure changes, application the energy principle the upstream conduit 
for any instant time yields** 


The velocity head for the upstream conduit included for convenience; 
also included are all form and surface-resistance losses between the reser- 
voir and point (1). Relationships similar Eq. (1) for other parts the surge 
system are given Appendix III. 

All acceleration heads are defined the conventional, one-dimensional 
formula 


(2) 


*The numbers parentheses indicate references listed Appendix 
**The symbols used are defined Appendix 
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The sign will depend upon the sign the derivative will 
appropriate use Eq. (2), Eq. (1) may written the form has 
ag con 
Similar relationships may developed for other Q’s. The diameter the 
main conduit has the same value for the lengths and L5. All corre- 
sponding dimensions the upstream and downstream risers and surge tanks 
are the same. 
The net head the turbine obtained from the relationship, 
Eq. (4), the comparatively small velocity-head lost the tail-race ig- 
nored. Moreover, assumed constant during any one solution. 
Application the continuity principle each tee requires that 
hea 
each instant, even when flow reversals take place. Differentiation Eq. 
(5) gives 
voi: 
similar expression may obtained differentiating Eq. (6). Equations 
the form indicated (7) are useful for eliminating the derivatives from 
equations the form indicated (3). Actual use shown Appendix III. 
The relationships between the water surfaces the surge tanks and the 
instantaneous flows out the tanks are obtained from the consideration (1) 
that inflow minus outflow equals the change storage during dt: 
(9) sin 
For cases when the water surface either tank drops below the level Eg, 
either Eq. (8) (9) replaced cor 
sys 
Head Losses 
All surface-resistance losses feet are evaluated the Darcy- 
Weisbach equation, 
ris 
(10) 
This equation applied not only the main conduit but also the risers. 
Surface-friction losses within the surge tanks are ignored because, ordinarily, 
they are negligibly small. 
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some installations part the total length the upstream conduit 
will lined with concrete and the remainder with steel. Provision for this 
has been made the machine program. Let length that lined with 
concrete and length that lined with steel, that Desig- 
nating and the friction factors for the concrete and steel lengths, re- 
spectively, and utilizing Eq. (10), the total surface-resistance losses from the 
reservoir point (1) are given 


Except for the tee losses, all form (or shape) losses feet are evaluat- 


(12) 


which assumed constant for particular flow-boundary configu- 
ration. Eq. (12) indicates that will have the same sign The total 


head loss* due various boundary forms from the reservoir point (1) 
given 


Strictly, all K’s Eq. (13) will have one set values when positive and 
will have different set values when negative (flow toward the reser- 
voir). this study the K’s are assigned one set values and are assumed 
remain constant during any one machine run. (If deemed necessary, the 
machine program can modified allow the K’s vary with the direction 
and also with the magnitude the velocity head.) 

Combining Eqs. (11) and (13), the total head loss from reservoir point 


(1) 


Head losses (other than tee losses) for other parts the main conduit are 
similarly determined, using analogous notation (Appendix III). 

evaluation the losses the risers and points (5) and (12), attention 
must given the direction and and whether not the water 
surfaces and are above below the elevation any instant under 


consideration. The total head loss for the upstream surge tank and riser 
system 


The first term Eq. (15) evaluates: (a) the loss velocity head when the 
riser delivering water the surge tank; (b) the velocity head the riser 
and the entrance loss when the surge tank delivering water the riser. 
The second term Eq. (15) evaluates the surface-friction loss the riser. 


*Strictly, hg, more than head loss, since includes the velocity head 
the upstream conduit. equal one all cases. 
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Eqs. and Appendix III indicate how and are evaluated. 
the two expressions for (Eq. allowance made whether 
not the total friction boundary has length something less than 

The third term Eq. (15) determines the head drop across the orifice and 
based upon the standard-orifice formula. this study taken being 
constant during any one machine run. Actually, even the orifice may 
considered thin-plate orifice such, likely that will vary 
changes owing the condition that the orifice plate probably not sufficiently 
distant from major flow disturbances most installations. However, may 
easily introduced into the machine program function appropriate 
variables, desired and data are available for realistic evaluation. Other- 
wise, may varied from one run another see what effect different 
values have the hydraulic transients, has been done the Oahe study. 
Head losses the downstream riser and surge tank are obtained the 
same manner used obtain 


Tee Losses 


determining the tee losses any instant time, the three Q’s each 
tee must known both magnitude and direction that instant. Evaluation 
the tee loss will then depend upon which the possible flow patterns indi- 
cated Fig. applies that instant. Provisions have been made the 
machine program for evaluating the tee losses each instant time for all 
flow patterns shown Fig. except the last two and j). (When flow 
patterns 2(i) and 2(j) occur, the tee loss assumed zero.) For these 
two cases, uncertain that appropriate data for evaluation exist. Should 
data available, the program could readily modified accommodate 
these two cases. Data are available(4) evaluate cases (a) through (h) 
Fig. Obviously, cases (e) through (h) are special cases (a) through (d). 

The tee losses feet are determined equations the general form, 


(16) 


the upstream tee the downstream tee) and upon 


whether the flow divergent convergent. The velocity head inserted 

stream tee) depending upon whether the main conduit flow directed toward 

(+) away (-) from the turbine. Eq. (16) replaced and for 


divergent flow and and for convergent flow. Four sets coefficients 


are required order obtain the piezometric-head changes each tee for 
all cases considered. (Table has been utilized for numerical evaluation 
when D,/D 2/3.) The coefficients include allowances for changes velocity 
head from conduit riser and vice versa. The details tee-loss determi- 
nation are given summarization Appendix III. 

With reference Fig. the tee losses are ignored, the piezometric 
heads points (1), (2), and (7) are assumed equal; that is, 
These three heads differ only when the tee losses are introduced. Thus, for 
example, 
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(17) 


(18) 


both these equations (applicable case (a) Fig. the sequence 


writing the subscripts for important from the standpoint avoiding con- 
fusion signs. Thus, indicates and not Pj. The same 


reasoning applies the downstream tee. 


all instances unsteady flow when the tee losses are taken equal 
are zero, the total energy heads the tees are not balanced because the 


instantaneous difference velocity heads between main conduit and riser 
(between points (1) and (2) Fig. for example). (The number times that 
this would occur during one total solution when the tee losses are evaluated 


all other instances during that same solution relatively small.) This un- 
balanced condition may assumed taken care orifice effects 


the riser, that is, the orifice coefficient. deemed necessary, special pro- 
visions can incorporated the machine program take care the 


velocity-head differences when this study, such provisions have not 
been made, because over-all results for the Oahe project, least, would not 


significantly affected. 


Gate-Time and Turbine-Performance Curves 


Transient effects within the hydraulic system are established turbine 
wicket-gate movements. order for the program work successfully, the 
gate opening per cent must known each instant time during the 
initial interval governor action, after which the gate opening assumed 


remain constant. This statement only partially true for surge-tank stabili- 
studies. (Discussion what involved stability studies presented 
another section the paper.) The gate-time curve given 


(19) 


time. Any single-valued function time may used, and may include 


dead and cushion time (for opening and closing gates, respectively.) Such 


curves are usually furnished the manufacturer the turbine. 
addition, turbine-performance curves are required. They enter the 
program table involving three variables, whose inter-relationship 


general mathematical form expressed 


H 


which function. For carrying out solutions for full-load rejection 


(20) 


maximum turbine design head and for full-load demand minimum 
turbine design head will usually necessary have tabular values 


and for values well above and well below Other- 


wise, the computer will stop because artificial bounds imposed lack 
tabular values regions computation which the analyst did not foresee 


overlooked. (What this means that when the computer reaches the end 
table and has instructions what next, stops.) 
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perform surge-tank stability studies along conventional lines, the 
following set turbine curves utilized instead the family curves indi- 
cated Eq. (20); 


which function and Bhp brake horsepower. This equation simply 
gives the relationship between and turbine discharge when the turbine 
power output constant and can include the effects varying efficiency. The 
use that has been made such curves and the difficulties experienced the 
computer study are discussed below. 


Machine Computation 


With the aim keeping machine-programming details practical mini- 
mum herein, Fig. has been especially prepared for the purposes this 
paper. The chart simply depicts the main elements problem solution but 
not the actual sequence machine operations.* show the actual sequence 
machine operations would require considerable detail, the inclusion 
which would tend obscure the main elements the problem. Fig. applies 
problem solutions that start from steady-state conditions, but the detailed 
machine program not limited. The equation numbers shown some 
the boxes refer the relationships given Appendix suggested that 
the reader first examine that appendix before proceeding with the remainder 
the paper. 

box (Fig. all lost heads are computed from the steady-state 
and data given boxes and and tailwater level adjusted 
the value must have for all heads balance exactly within the machine for 
the given reservoir elevation (box 2). Alternatively, the tailwater level 
may kept constant, and may then adjusted for balancing steady- 
state heads. This balancing process insures that the computer starts off with 
totally consistent set values for the variables involved. 

Once steady-state conditions are balanced within the machine, transient 
set applying the initial time interval (box 5), which fixes the new 
gate setting from box The new turn fixes the new from 
the turbine-performance curves (box 10), using the value the begin- 
ning the time interval. result, all other variables are changed, and 
the machine will compute their values each step the transient develops. 
The initial arbitrarily selected the analyst and can small 
seems warranted. Values low 0.01 second have been used start the 
transient; most cases, initial 0.10 has been used. The computer can 
instructed through control information (box use values 0.1 
total time 1.00 and use time intervals one second thereafter 
to, say, seconds, after which would use 5-second intervals until the total 
time for the machine run (given box reached. These figures are mere- 
illustrative. Any reasonable set time intervals may used, including 
possible use constant the end the run. However, the method 
solving differential equations (discussed below) will not work expeditiously 
the assigned time intervals for step solutions are too large. 


*Complete details, including all flow charts and the machine codes, are given 
the machine program itself, which being prepared separate 
document for distribution limited basis. 
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Read in all physical data: L., bs, 
E,. Eg. and K, for Surge Tanks. 
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Read in main program instructions and al! tables: 
1, =F, (Mg); Tables for Fy. fa. fg, If different; 


Table for Tee Losges; H=F(G,Q,); G=F,(t) 


Read in contro! information: Tee losses computed or ignored? Riser acceleration 
computed or ignored? Initial time intervalSt; Total time for run; Initial valves 
of Qs. Vatues of at beginning and end of governor action. 
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With reference Fig. computations are carried out the loop, boxes 
through and back (but not precisely the sequence shown), with 
out results box occurring whenever all the conditions for solution for 
each time step are satisfied. The computer may work around this loop the 
times during any time step before print-out results occurs. The computer 
may print out all any number the quantities indicated box 26. box 
the indicated relationship is, course, never true. This simply 
trick which used save some programming steps for rejection cases 

erv 
the computer determines the head losses (box 13) and tee losses (box 15) 
every step. manual methods, this would entail tremendous amount 
work. calculate all the each step requires the determination 
five Reynolds Numbers and linear interpolations the tables for six 
factors (box 4). Tee losses are evaluated accordance with all the con- 
ditions specified Appendix III and linear interpolation table having 
four sets coefficients. This degree refinement surge-problem compu 
tations practicable because the machine costs involved are comparatively 
con 

The choice whether not include riser acceleration (box 16) not that 
necessary iuture applications. The option was incorporated this 
program order give direct indication what effects this item would 
have final results. option computation tee losses has also been 
incorporated the machine program, but the option not indicated Fig. 
(box 15). Here again, future applications, there would necessity for 
ignoring tee losses. long the form the tee-loss equations and the 
method evaluation not change, any set numerical values may 
ized for the tee-loss coefficients the machine program. The 


nated 


Machine Solution Differential Equations 


Knowing the values the primary variables Q1, Q3, Q5, Ws1, and 
the beginning time step, the derivatives these five vari- 
ables (boxes 14, 22, and Fig. may readily determined for time 
ty. and are not designated primary variables, since their values 
are easily obtained algebraic subtractions through use continuity Eqs. 
(5) and (6)). The values the primary variables the end the time step, 
are obtained process high-order approximations. (10) 

For ease writing let equal any one the five primary variables for 
the purposes this section the paper. Then the values any one 
quantities the beginning, mid-point, and end the time step are origir 
nated yj, and yg, respectively. Roman numerals are utilized 
nate approximate values and Thus, indicates the first proce 


known. first approximation made utilizing the value the 


just 


Having first approximation the five primary variables, indicated *The 
Eq. (22), new values (temporary) their five derivatives are then computed| 
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exactly the same way was done time these new derivatives are 
nated The next step yields better approximation making use 
for 
any the the two derivatives now hand, that 
uter 
Having set given Eq. (23) aside, the computer cuts the original time 


interval half and determines exactly the same manner was 
calculated, with the result that 


q 


vely Still another approximation for obtained applying the same method 

computation the second half the time interval (from tg), 

not that 
Fig. 
util- 

The final values all primary variables time are taken as* 

ari- 


time But before these final values are assigned yo, relative error tests are 
lues made. The relative error defined 


Eqs. 

these pre-assigned margin relative error, the computer cuts the 


original time interval half and goes through the whole procedure (em- 
bodied Eqs. (22) through (29)) all over again. The machine will repeat this 
process until (For the Oahe project, Was given value 
prevent the computer from taking unnecessarily small time increments, 
the machine instructed double the original time interval whenever 
Ry| another pre-assigned margin relative error. (For the Oahe 

project was assigned value that the procedure 

just outlined, solutions are advanced time step time step and automatic 
error control achieved each step. 


*The factor three Eq. (28) based upon theoretical considerations. (10) 
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those made independently engineers Sverdrup and Parcel, Inc., firm 


The necessity for such high-order approximations for solutions differ- 
ential equations has not been investigated, because the machine time required 
per total solution reasonable and because office time could not spared 
for the mathematical explorations that would required. believed that 
the machine method solving differential equations for surge problems, 
discussed herein, can also applied profitably machine solutions 
variety other engineering problems. 


Application Machine Methods Oahe Surge System 


Conditions and Data for Numerical Solutions 


Oahe Dam, located about six miles upstream Pierre, South Dakota, 
key project the system multiple-purpose reservoirs the Missouri 
River.(5,6) Each the scheduled seven generator units Oahe will have 
nameplate capacity 85,000 kilowatts, and hydro-power will delivered 
the generators Francis-type turbines. The hydraulic surge system 
ing any one the seven units schematically shown Fig. The 
the upstream conduit different for each unit due physical layout. 
Machine results have been obtained only for unit number for which 
3,820 feet all cases, except cases and Tables and These| 
two tables and Table give the initial conditions (steady state) and other 
pertinent data upon which manual and all machine computations for full-load 
rejection, full-load demand, and surge-tank stability analysis are based. 
References manual computations all tables and figures herewith apply 


under contract with the Omaha District, Corps Engineers, for design 
Oahe power plant and surge tanks. 

All computations are based upon: surge tank diameters feet with 
interconnection between tanks; riser diameters feet; main conduit and 
penstock diameter feet; riser height feet; orifice diameters 
12.33 feet; distance between risers feet (except for cases and 
Ile); and other data given Tables and With feet, the surge 
tanks will, course, not fit concentrically over the risers. Hence, for 
and the distance was doubled order permit concentric 
rangement surge tanks and risers, desired. From hydraulic stability 
standpoint, the distance between tanks should probably small 
cable. The eccentric arrangement with feet has been adopted for 
Oahe project. 

All friction factors for both manual and machine computations have been 
evaluated from data given Tables and When friction 
losses for the smallest diameter (16 feet) are negligibly small, and hence, 
evaluate all tee losses for computer solutions. Manual computations all 
cases are based upon ignoring acceleration effects the risers 
length but tee losses have been evaluated using Vogel’s will 
noted from Tables and that for manual computations was 3,580 
instead 3,820 feet and that few other items for these computations differ 
slightly from those used for machine solutions. The reason for this that 
computer runs were made much later data, after some Oahe design 


changes had been effected. Other conditions and data well results are 
discussed below. 
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TABLE 5 


Tee Loss Coefficients* 


surge 
cases Discharge Ratio 
bility j Sconduit 
lor the 
0.2 0.85 | -0.87 0.9 0.68 0.83 0.15 
n all 0.8 $.20 -0.75 5.95 1.0 2.17 | 1.0 
3 1.0 7.0 | -0.65 8.50 | 0.70 2.45 | 1.75 
the 
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Results for Full-Load Rejection 


With reference Table all cases are based upon rejecting full load 
the turbine, 142,000 horsepower output, maximum net head the turbine 
203 feet for steady state. For this condition the gate opening restricted 
66.3 per cent for all computer solutions (70 per cent for manual compu- 
tations) order avoid overload the generator, and closure the gates 
achieved 3.31 seconds (computer) uniform rate per cent per 
Turbine-performance values (Eq. (20)) have been obtained linear 
table the form indicated Table The actual turbine- 
performance table that was used machine computations gave values for 
and 240 feet). 
Figs. and indicate some the results that were obtained for the re- 
jection cases. Fig. depicts the water surface elevation the downstream 
surge tank and the difference water surface elevation between the two tanks 
function time for all cases indicated Table Case and Manual 
Computations are based upon the same data and conditions, except that for the 


manual computations over 200 feet shorter, initial and are slight- 
higher, and acceleration heads the risers and the length are not 
computed. The basis computation for all other cases shown Figs. and 
the same for Case except that for: 


Case Ib, acceleration heads the risers have not been computed. 
Case Ic, the orifice coefficient 0.8 instead 0.7. 
Case Id, the orifice coefficient 0.6 instead 0.7. 
Case Ie, the length feet instead feet, and being some- 
what smaller consequence. 


Examination the curves Fig. indicates successful application 
machine methods rejection problems. noted that: (a) the numerical 
value the orifice coefficient has the most pronounced effect upon the water 
surface elevations the tanks and (b) comparison with the results manu- 
computations, the difference water surface elevations between tanks 
substantially greater for all computer cases during the first quarter cycle. 
Actual machine-computation time for each the cases shown Fig. aver- 
aged slightly over one hour. The values for Case Fig. are more 
less typical those for all other computer cases indicated Fig. 


Results for Full-Load Demand 


With reference Table all cases are based upon picking load when 
the reservoir elevation minimum design level (Ep 1540). The rate 

the wicket gates linear all cases, the rate being per cent per 

second (governor time seconds) for all except Ila. For Case Ila, the 

time seconds. Applicability Table load-demand cases 

the same for the load-rejection cases previously discussed. Case and 

Manual Computations (Table are based upon the same data and conditions 

except that for the latter: over 200 feet shorter; initial steady-state 

zero (instead 5.0); gate opening end governor action per cent 

100 per cent); tailwater elevation 1420 (instead 1421.7); and 

heads the risers and the length are not All 

cases are based upon initial steady-state discharge 390 

(corresponding initial 5.0) instead zero. comparatively small 
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‘cate oscillations rather large amplitude between tanks for short period, 


finite value the start seems necessary because the machine 
methods solving differential equations that has been utilized. smaller 
value than 390 c.f.s. might have been used successfully. However, provided 
that the initial small comparison the full-load the over-all re- 
sults are not likely significantly affected. 

Figs. and indicate some the results that were obtained for the 
demand cases. data shown under Case Table also apply Cases 
and except that the surge tank base elevation 1515 does not ap- 
ply Under Case the surge tank water surface drops down into 
the riser, but under Case II-2, sufficiently lower (about feet) 
the water surface does not recede into the riser. (With this small difference, 
the change ignored.) The basis computation for all other computer 


cases the same for Case II-2 except that for: 


DIFF FT. 


Case Ila, the governor time seconds instead seconds, and 
with 1515 the water surface drops down into the riser. 

Case acceleration heads the risers have not been computed. 

Case Ild, the orifice coefficient 0.6 instead 0.7. 

Case Ile, feet instead feet thus causing and 
somewhat smaller, and 1.00 (re-entry tube) instead 0.5. 


machine methods have been applied successfully also full-load demand 
problems. The question whether not permit the water surface drop 
down into the riser actual operation immaterial insofar the purpose 
this paper concerned. Results are now available basis for decision 
that respect. will noted from Figs. and that machine results indi- 


Comparison manual and machine results Fig. demonstrates that 


when the water surface permitted recede into the risers. Cases II-1 and 
are directly comparable except for governor time, and hence results for 
these two cases are plotted separately for comparison Figs. and 
Doubling the governor time seems make little difference over-all 
results for the Oahe surge tanks such; that is, low point water surface 
elevation the risers (Fig. about the same and the timing therefor 
also about the same. But, course, with slower governor, pickup the 
turbine load materially delayed, comparison values Fig. indi- 
cates. interesting note from this figure that except for slight lag 
time, the results Case Ila are almost carbon copy those Case II-1. 
This demonstrates that the computer does the job with slavish adherence 
instructions, and, consequence, the analyst cannot afford slight 
overlook any pertinent detail the preparation the machine program. 
machine-computation time for each the demand cases discussed above 
averaged about one and one-half hours. 


Surge-Tank Stability Analyses 


normal practice, stability analyses are made for small power-demand 
changes minimum design head the turbine the region drooping ef- 
ficiency (between gate opening for maximum efficiency and full gate.) Ideally, 
what desired directly from the initial steady-state power 
ceed directly from the point 75,000, 120) point 
78,000, 119.6), the new steady state for constant turbine output 
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78,000 horsepower. However, this not possible, for sudden gate 
movements response rapid governor action set transients the 
surge system. These transients will die out reasonable length time, 
the system properly designed. conventional methods, assumed for 
purposes analysis that once entry onto the curve made, 
point that may be, the power delivered the turbine will have constant 
value thereafter. The oscillations andG which occur 
constant are then determined manual methods arithmetic inte- 
gration. Such methods are necessarily based upon simplifying 
result which the total transient problem treated one 
multiple parts. essence, the hydraulic system divorced from the inter- 
connected mechanical and electrical system for the practical purpose ana- 
lyzing the hydraulic transients such. Use this same approach 
lutions computer (employing the method solving differential equations 
previously discussed) has been unsuccessful, but only the sense that the 
computer has failed produce the same results those obtained manual 
computations. another sense the computer study has been quite fruitful 
that has forced reappraisal the total “transient” problem and has 
lighted the shortcomings manual methods analysis. Discussion the 
need for comprehensive analysis power plant transients integrable 


whole presented another section the paper. 

With reference Table all computer runs were made starting with 
Bhp, 75,000 for initial steady state net head the turbine either 
130 feet (the former value being the minimum design head for the Oahe 
turbine.) The gate-time curves indicated Table are shown Figs. 


10, which the initial steady-state gate opening (line one Table 8). 


Curve No. Fig. represents part the full-range gate-time curve 
furnished the Oahe turbine manufacturer. Its relation Curves No, 
and indicated Fig. 10. Curves and are entirely arbitrary. 
for all computer runs are exactly the same all respects except for 
initial values (line Table 8), gate curves, and indicated inclusion 
exclusion acceleration heads the risers (line Table 8). 
The first computer run that was made for stability analysis surge 
has been designated Test No. and involved use Gate-Time Curve 
(Table 8). Data for this test are exactly the same for Manual Compu- 
tations, except for the inclusion acceleration heads and the 
The path taken the computer trying reach 78,000 
indicated Test No. Fig. 11. Full gate was reached 1.1 seconds be- 
fore was reached. following the paths indicated Tests No. 
and Fig. 11, the computer was instructed follow the gate-time curve 
given (in tabular form) until reached the curve Bhpg, after which con- 
trol for development the transient would revert that curve. (That is, 
thereafter the function represented Eq. (21) instead Eq. (20) would have 
satisfied well other relationships.) Also, the computer was pro- 
grammed stop whenever reached full gate, since oscillations past full 
gate are unacceptable. 
Test No. having failed, Test No. was started from region 130) 
where the efficiency curve flat and turbine efficiency practical constant 
down 120. Again, the computer followed devious path, reaching 
1.36 seconds (Fig. 11), after which continued along the curve Bhpg full 
gate. With the thought that, perhaps, riser acceleration effects caused the 
devious path for Test No. Test No. was started from the same point 
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130) with the effects riser acceleration excluded. This time, the path 

followed the computer was satisfactory until the Bhpg curve was reached 
0.75 seconds (Fig. 11), but after that the computer followed the curve 

upward full gate once again. These failures establish transient (simi- 
lar that obtained Manual Computations) that would produce oscillations 
and along the constant Bhp curve 78,000 without any time 
reaching full gate not indicate that the Oahe surge system unstable 
the region investigated. Rather, they indicate that any transient set the 
manner conventional methods analysis likely resemble the actual 
transient only accident, for these methods ignore (with justification*) the 
oscillations the gates with time due other effects the turbine and 
governor end the surge system. anticipation any doubts that may arise 
concerning the adequacy the machine program, sufficient data are given 
this paper for the reader assure himself that the computer has done the job 
accordance with instructions. Machine print-out results for Test No. 
are given the upper part Table These data together with the equations 
Appendix III, data given Tables and and gate-time Curve No. 
Fig. will permit checks machine arithmetic relatively short time. 
will noted that the machine does not print out values and they 
may obtained separate computation. 

When load change called for external electrical power require- 
ments, oscillations are set the turbine governor system itself during the 
first several seconds. virtue interconnection, the gates also oscillate. 
effort represent least the first one-half cycle these oscillations, 
gate-time Curves No. and (Fig. 10) were invented, using the “dead” 
and “cushion” ends the manufacturer’s curve (of which Curve No. 
part) form the peak regions and the tail-ends. Application these gate- 
time curves then gave the results shown Fig. 12. all three cases the 
computer started and stayed within region which turbine efficiency 
practically constant. following the path indicated Curve No. Fig. 12, 
the computer reached 0.8 second, after which proceeded along 
full gate the same manner before. (Machine print-out some the 
results based use Curve No. given the lower half Table 9). 
this point the study program change was made which the machine was 
instructed ignore the first intersection with and follow gate-time 
Curve No. (Fig. 10) until reached second time, after which control 
would revert For this condition the machine stopped 1.2 seconds 
(near the end gate-time Curve No. shown Fig. 10), and this test was 
abandoned. This same procedure was applied gate-time Curves No. and 
with the results shown Fig. 12. following the path marked Curve No. 
the computer followed instructions perfectly, reaching second time 
0.925 second, after which followed downward the right and stopped 
1.125 seconds. stopped this point because had reached the end 
the table giving function along the curve Bhpg, the end the table 
having been purposely set limit 140. Gate-time Curve No. was 
applied with the aim having the path loop back and intersect the 
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—— 


*Justification exists because, normally, these oscillations cannot described 
adequately without bringing into the problem all the relationships that tie the 
hydraulic, mechanical and electrical systems together. Manual analyses 
this basis would totally impracticable. 


’ 
| 
3 
) 
© 
Ae 
_ 4 
t 
3 
J 
— 
— 
— 
+ 
NJ 
— 
— 
— 


POWER PLANT SURGE TANK ANALYSIS STABILITY STUDY TURBINE UNIT NO. 


TABLE 


Tailwater Elevation 1420.00 


Gate Curve No. 


April, 1959 
=A 


1553. 


oir Elevation 


N 


ion 


Tailwater Elevat 


130 


124 
124 
123 
123 


ean 


42/5810 


LEGEND 


q 


q | 
a 


0099 0099 
<= 


a 
a 
| 
te 
re 
3 
= 


opposite direction from that which had been indicated Curve No. The 
computer started return loop along Curve No. indicated the points resul 
0.9 and 1.0 second Fig. and then stopped 1.0 second (the end 
gate-time Curve No. Fig. 10). this point was decided that there was 
altogether too much arbitrary judgment being exercised, and further comput- 

solut: 


studies stability were not undertaken. 


Machine-computation time for all these tests was comparatively small. Thet 
However, good deal time was spent analyzing machine results 
trying establish transient using function the form given Eq. (21) turba 
concluded from the experience thus gained that analyses Obvic 
utilizing this procedure yield artificial results and that, consequence, 
the results conventional methods analysis surge-tank stability 


the lines previously described are incomplete and inconclusive. This does 


not raise imply any reflection the excellent work that has been done 
connection with hydraulic design power plants. For, has long been 
nized that simplifications are practical necessity order render power the 
plant “transient” problems tractable for solution manual methods 
reasonable length time. However, successful adaptation for computer so- numb 
lutions very complex problems many fields engineering and science histo 
now raises the question the necessity for continuing use simplified 
procedures that yield results that are not entirely satisfactory. 
mach 

Some Ideas for Comprehensive Study Digital Computer the 

The ideas presented this section the paper have been extracted from Aside 

unpublished “Proposal for Comprehensive Study Digital Computer powe 
Hydraulic and Governing Transients Affecting Design and Operation 
Power Plant System”, prepared the writer his own initiative and all 
mitted the Corps Engineers for consideration. This proposal includes meth 
the derivation necessary equations and suggested methods machine nary 
sis and computation for undertaking comprehensive study, which only ing 
briefly discussed below because space limitations. The only reason for in- the 
viting attention the proposed study this time indicate the engi- withi 
neering profession that efforts are being made overcome the limitations write 


conventional methods. scope 
Formidable limitations time and means have required now that 
power plant system treated divisible whole for purposes analysis 
hydraulic and governing transients induced external power-load vari- 
ations. (As used herein, the term “power plant system” includes the surge 
system previously defined and one power unit, turbine and generator, sup- 
plying load independent basis.) With the advent powerful, high-speed 
computers, appears that, for the first time, practicable analyze the 
total “transient” problem “feedback” system interwoven effects rather 
than problem having separable, multiple parts. Studies water hammer, 
surge-tank effects, turbine-speed regulation, and governing stability need 
longer made essentially separate basis computer program can 
developed handle the total problem such. Recent analog-computer 
studies Koenig and Knudtson(9) are only step toward the ultimate 
for these studies involve simplifications which tend cast some doubt the 
validity results obtained thereby. 
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appears that studies the “total problem” can made with conclusive 
results only when the relationships controlling hydraulic behavior the 
conduit, penstock, and surge-tank system are solved simultaneously with the 
connecting relationships between turbine output, the governor, machine in- 
ertia and speed, and external power load. requires the simultaneous 
solution great number equations, short time step short time step. 
The time step would probably have the order hundredths 
second (or less some cases) order account for all pressure dis- 
turbances and reflections* generated during the existence transient. 
Obviously, such approach totally impractical, not impossible, for 
manual computation, but, apparently, not for modern digital computer hav- 
ing the capability and computational speed comparable with those IBM 
704 computer. 

has been roughly estimated the writer that the machine time required 
per total solution such computer would well within practical range, 
the necessary machine program can successfully developed. total so- 
lution would yield time histories velocity and piezometric heads large 
number pre-determined points within the surge system well time 
histories turbine power output, turbine discharge and head, turbine and 
generator speed for any assumed given external power-load variation with 
time applied the power unit under study. Important quantities such 
surge tank size, orifice characteristics, penstock and conduit characteristics, 
machine inertia, governor characteristics, etc., would all introduced 
the total problem parameters. the proposed methods, water-hammer 
and mass-oscillation phenomena would not treated separate entities. 
Aside from the basic assumption essentially independent operation the 
power unit, other simplifications important significance appear 
necessary for the methods solution that have been proposed. accomplish 
all this digital computer requires the simultaneous solution numerical 
methods partial-differential equations within the surge system and ordi- 
nary differential equations the turbine end the system. Methods for do- 
ing this are apparently available. remains determined whether not 
the fundamental partial-differential equations describing hydraulic behavior 
within the surge system can solved computer practical basis. The 
writer believes that this can done, but machine programming job great 
scope will have performed before success achieved. 


CONCLUSIONS 


machine program for the IBM 650 computer now available for solving 
load-rejection and full-load demand problems surge system comparable 
that the Oahe project one two hours machine-computation time 
per case. The machine time per solution for installations having only one 
riser connected the main conduit will materially less. The program de- 
scribed herein may easily adapted for such cases. The relatively small 
time required per machine solution and the comparatively small costs in- 
volved make feasible enlarge the scope design studies surge 
systems well beyond the limits ordinarily imposed manual methods so- 
examining the results number machine solutions readily 


*In the study that envisioned, effects the compressibility water and the 
elasticity flow boundaries would included. 
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obtained result changing one two major factors from one machine 
run another, better understanding problem phenomena for particular 
installation will obtained. consequence, construction savings may 
achieved some cases. 

digital computer program which will yield solutions for transient effects 
resulting from small-load changes from one steady-state power output toan- 
other not yet available despite considerable efforts develop such 
program. may entirely possible develop such program utilizing 


methods solution other than those described herein. However, appears 
that any method which ignores the applicable relationships between turbine, 
governor, machine speed and inertia, and external power load for the purpose 


analyzing the hydraulic transients such bound yield inconclusive 
results. The premise that such simplifications give conservative results in- 


sofar surge-tank stability concerned may true but the degree con- 
servatism not known application particular cases. For long time 
need has existed for analyzing the total “transient” problem without the major 
simplifications that are ordinarily made. When only manual methods 
sis and computation were available, fulfillment this need practical 
basis was not possible. With the advent powerful, high-speed computers 
the prospects are now excellent that this need can wholly satisfied. Some 
efforts have already been made toward the ultimate goal adapting the total 
“transient” problem power plant for solution digital computer such 
way that machine results for water hammer, surge-tank effects, turbine- 
speed regulation, and governing stability would delivered one package 
few hours machine-computation time. successful adaptation ulti- 
mately achieved, many benefits and possible construction savings will 
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Notation 
Cross-sectional area sq. ft. conduit and penstock. 
Area orifice sq. ft. 
Area riser sq. ft. 


4 
hp 
| 
¥ 


Area surge tank sq. ft. 
Bhp Brake horsepower (output turbine) 
Coefficient for risers, 


Cp; Coefficients loss for diverging flow tees, function dis- 
charge ratio (Qriser/Qconduit) 


Coefficients loss for converging flow tees, function dis- 
charge ratio (Qriser/Qconduit) 


Diameter feet conduit and penstock. 

Diameter surge tank feet. 

Absolute roughness length feet for concrete. 

Absolute roughness length feet for steel. 

Pre-assigned margin relative error. 

Elevation reservoir feet above M.S.L. 

Elevation base surge tanks feet above M.S.L. 

Elevation centerline turbine distributor. 

Darcy-Weisbach friction factor, function Reynolds number. 

Friction factors upstream conduit. 

fo, Friction factors upstream and downstream risers, respective- 
ly. 

Friction factor penstock between risers (Length Lg). 

Friction factor penstock (Length Ls). 


Functions Reynolds Number for concrete and steel. 


Gate opening per cent for turbine wicket gates. 

Initial gate opening for steady state. 

Surface resistance loss feet. 

Form loss feet. 

loss feet (may include all types losses and velocity 
ead). 
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Head loss feet upstream conduit length (includes 
all form and surface-resistance losses well velocity 

head), defined Eq. 


Head losses upstream and downstream risers, 
defined Eqs. and 


Acceleration heads upstream and downstream risers, re- 
spectively. 
Acceleration head penstock length Ls. 
constant, 


Constants which apply upstream and downstream risers, re- 
spectively, defined Eqs. and 


Entrance loss coefficient for upstream conduit. 
Entrance loss coefficient, surge tank riser. 
Loss coefficient for trashracks. 
Loss coefficient for all gate slots. 
Kts Loss coefficient for all transitions. 
Loss coefficient for all horizontal bends. 
Loss coefficient for all vertical bends. 
Kot Coefficient account for all other form losses. 
Lengths feet main conduit (see Fig. 1). 
Portion lined with concrete. 


Tj 


Piezometric head feet above M.S.L. 
Piezometric heads indicated points (see Fig. 1). 
Discharge c.f.s. upstream conduit. 
Discharges c.f.s. upstream and downstream 
risers, respectively. 
Discharge c.f.s. penstock length (Fig. 1). 
Discharge c.f.s. penstock length (Fig. 1). 
Indicates absolute value (ignore sign). 
Height riser feet (Fig. 1). 
Time interval. 
Tee loss feet. 
Tee loss feet between points and (Fig. 1). 
Tee losses feet between indicated points. 
Tailwater elevation feet above M.S.L. 
Vo, Velocities ft./sec. upstream and downstream 
risers, respectively. 
Velocity ft./sec. length Ls. 
Used represent any one variables, Qs, 


and 


Water surface elevation (M.S.L.) upstream surge 
tank. 

Water surface elevation (M.S.L.) downstream surge 
tank. 

Kinematic viscosity water feet squared per 
second. 
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APPENDIX 
Mathematical Relationships 
Notation 
Detailed notation given Appendix general, the following notation 
utilized, based upon the pound-foot-second system: 
Coefficient (riser, orifice, Coefficient for losses 
Absolute roughness length Piezometric head (M.S.L.) 


ay 
jal 
= 


| 


hy, Lost head Velocity 
Velocity head Water surface elevation 


Sign Convention 


The net head the turbine and all piezometric heads are always positive. 
The signs all lost heads and acceleration heads are automatically assigned 
the equations that define them. The signs all velocities will the same 
those the corresponding discharges, since AV. The sign convention 
for all discharges follows: 


All flows directed toward the turbine are given positive (+) sign. 


All flows directed away from the turbine are given negative (-) sign. 


Flows the risers are positive when directed downward and negative 
when directed upward into the surge tanks. 


Mathematical Relationships 


The physical equations and functions required solve surge problems for 
the system shown Fig. are given this appendix. The major variables 
which vary with time tare: all Q’s, all V’s, all P’s, all f’s, 
Cp, Cp, Ce, and Such quantities Co, etc., may enter the 
machine program parameters; that is, such quantities are constant during 
any one run, but may varied from one run another. 
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Substituting from and Eq. becomes 


Using and eliminate Pg, from obtain 


+ 5 
Solving and simultaneously yields 


evaluate first necessary determine the head losses and 


Head Losses are determined Eqs. through follows: 


4 
a 
= 3 
~ 
“4 
i 


Tables values for and are needed for machine runs; see, for ex- 
ample, Tables and herewith. Similarly, for all other friction factors, 
fo, all which are functions and same tables may used. 
+ 
13) 
15) 
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A 


long the form the equations and the method evaluation tee 
losses not change, any set numerical values may used for the 
coefficients the machine program. 


Tee Losses are evaluated accordance with the following summarization 
Cases through (Eqs. through The coefficients Cp, 
and are evaluated from tables; see, for example, Table herewith. 
For ease writing, let hy; hy3 and and recall that 
Case Tee losses ignored. 
Case and are both positive and negative. 
and are both positive and negative. 
and are negative and negative. 
age D D 8,9 D D Vs (11-34) a(t 
hy = -C hy (mI 35) 
{ 
and are positive and positive. 
1,7 8,14 (I-36) | 
: 
( 
and are negative and positive. 
FY. 
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SURGE TANKS 


Acceleration heads are defined follows: 


h = 


when 


s 


Water surfaces surge tanks depend upon following equations well 


others: 


A_, when 


= is 


when 


when 


Turbine end the surge system; the following relationships must also 


satisfied: 


Gives the gate opening function time. Either given manufacturer 
must The function enters the machine program table; 


hence, may have any form. 


Turbine-performance curves enter the machine program table; see 


for example, Table herewith. program, and are the independent 


ion 
hat 
ble 
{ 


variables and the dependent variable; that is, computed when and 
are known. 


Bhp constant 


This function used only for stability studies; see, for example, Figs. and 
herewith. These curves enter the program tables. 
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TWO METHODS COMPUTE WATER SURFACE PROFILES 


SYNOPSIS 


Two methods derived the authors their work with the Bureau 
Reclamation are presented for the computation water surface profiles used 
for the development tail water rating curves and backwater profiles. Ex- 
amples each method are included which show the computational procedure. 


INTRODUCTION 


The derivation water surface profiles serves two major purposes the 
Bureau Reclamation; one—to determine tail water rating curves, and two— 
trace the backwater curve above adam. Tail water curves are used the 

power-plants, pumping plants, and energy dissipators such still- 

ing basins. the design large dams, these curves also furnish useful 

information for making stability and stress analyses. primary use the 

backwater curve above dam for reservoir land acquisitions and 

Backwater information also necessary the design new fa- 
cilities such pumping plants, bridges, powerplants, other structures 
located along above the reservoir. 
The methods for computing water surface profiles developed herein are 
referred Method and Method Method can used computing 

profiles for channels which the flow distance between sections equal for 

segments the flow. Method differs from Method that es- 

pecially adaptable conditions where the flow travel paths are longer for the 


proper hydraulic analyses the conveyance capacities each cross- 
sectional subdivision. Eddy loss corrections and velocity head changes can 
also taken into account either method. 


Note: Discussion open until September 1959. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 1997 
part the copyrighted Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. April, 1959. 

Hydr. Engr., Bureau Reclamation, Denver, Colo. 
Hydr. Engr., Bureau Reclamation, Denver, Colo. 


main channel than for the overbank portion. Both methods, however, are based 
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Theory Other 
The underlying theory each method based Bernoulli’s Energy The 
evaluate 


Equation which repeated for convenience the following discussion its 
application water surface profile computations. 
The conditions open channel reach steady, non-uniform flow are 
diagrammatically illustrated Fig. The energy balance equation and 

ten: 


where where 
streambed elevation referenced given datum (downstream 
section) 
hy, velocity head downstream section 
Sip 
eddy 
The 
terms 
The tota 
equation 
where 
using 
sir 
ENERGY OPEN-CHANNEL FLOW 
methods 
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WATER SURFACE PROFILES 


Other losses eddy, bend, and bridge pier losses. 


The two methods use basically the following hydraulic relationships 
evaluate the various components the energy equation. 


2 


and 


where 

discharge cubic feet per second 

cross-sectional area square feet 

cross-sectional velocity feet per second 

gravitational acceleration feet per second squared 


reach length between sections feet 
roughness coefficient representing the reach 


hydraulic radius equal cross-sectional area divided the wetted 
perimeter 


noted from above that Manning’s equation used compute 
Generally, Manning’s discharge equation written follows: 


2 a 


where the slope the energy gradient. 


The friction head can determined applying Manning’s equation 


hg, eddy loss, bend loss, and bridge pier loss. 
terms the friction slope, 


2 
The total friction head for reach then can resolved the following 


equation: 


where and Sfo are the friction slopes computed the preceding equation 
the representative hydraulic properties: and for Sections and 


head and vice versa open channel flow. While there are, present, 


The total losses the energy balance equation consist the friction head, 


simplified and practical concept eddy loss can best understood 
Studying the energy losses incident the conversion kinetic head static 


methods for measuring these losses, known they occur and manifest 
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themselves usually the form turbulence. Observations channel flow 
where there are abrupt changes velocity clearly display turbulence the 
presence boils and vortexes. This particularly evident when there Method 
significant change cross-sectional area. This 

has been general practice the Bureau Reclamation compute the tations. 
eddy loss 10-per cent correction the difference velocity heads (Ahy) the 
when the static head converted kinetic head and per cent for cable 
condition. Thus, when the value hy, positive, 10-per cent cor- channel 
rection applied, and for negative values 50-per cent correction used. locity 

Departures from the 50-per cent value quoted above may necessary, rection: 
particularly, cases where extreme changes velocity occur the channel represe 
reach. This condition results transitions from very narrow section the mai 


wide section. Large standing waves are formed, which indicates the flow Refe 
bordering critical stage. the flow progresses downstream, the energy termine 
dissipated turbulence and upon reaching the downstream section, reverts know 


more tranquil state. some these cases, the eddy loss could reason- 
ably increased much 100 per cent the change velocity the floo 
always applying 50-per cent correction the computations, situation those 
can result showing higher water surface elevation the downstream section 


dg, adverse water surface slope) which seldom en- The 
countered natural open channel flow. correct this condition, ad- and 
justment can made increasing the percentage applied. Eddy losses, equal 
best, are only estimates based the judgment the hydraulic engineer. channel 

search the literature reveals that there have been formulas de- Vs. 
veloped yet for thorough determination bend losses. number formu- 


las resulting from various studies are listed following page. 


channel width 


0.21 inner radius 


from Yarnell and Woodward(1) 


“n” increased 0.001 for each 20° curvature suggested Scobey(2) 


0.0256 width stream 


then 
graphic radius curvature 


Boer and Urick(4) 


There are set criteria which have been established make practical 
evaluation bend losses. Two methods are suggested for the determination 
these losses; first, they could computed from one the formulas listed 
above; and second, “n” values could modified include bend loss effects. 
Ordinarily, such losses determined the first method are very theoretical 
and require analyses usually beyond the limits practical dimensions the The 
collected field data. recommended, therefore, that “n” coefficients the 
assigned values include any significant effects bend losses. The 
The hydraulic conditions bridge crossings require careful analysis 
determine bridge pier losses. King’s Handbook(5) recommends contracted 
and enlargement coefficients 0.5 and 1.0, respectively, applied the 
velocity heads the approach, bridge, and emergent section. Nagler, 
and Johnson(9) have developed formulas for computing the head 
loss through bridges which may applied under certain conditions. 
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Methods 


Method 


This trial and error procedure which involves step-by-step compu- 
tations. The method has very wide scope application and can used 
the general case involving steady, non-uniform flow. particularly appli- 
cable the irregular channel which the cross section consists main 
channel and separate overbank areas having individual “n” coefficients. Ve- 
locity head changes are taken into account weighting process and cor- 
rections can included for eddy losses within the reach. Reach lengths 
representing the flow path between sections, however, are assumed equal for 
the main channel and overbank areas. 

Referring the diagrammatic sketch Fig. the objective de- 
termine the water surface elevation for given discharge Section when 


When the discharge exceeds main channel capacity, the flow also inundates 
the flood plain overbank area. channel bed conditions are different from 
those the overbank areas, separate “n” values are assigned, and the cross 


section subdivided for proper hydraulic analysis. 

The areas and hydraulic radii are computed separately for each subdivision 
and tabulated shown for Section Table The conveyance capacity, Kg, 
equal 1.486/n AR2 also shown this table, computed for the main 
channel and overbank using the corresponding “n” value. Conveyance curves, 
vs. elevation, are plotted for each section; example shown Fig. 

this method, the friction slope, determined the following 
analysis: 


2 a 
since 
1.486 
then 
solving 
ted 
the The value used the above equation the total conveyance capacities 
the main channel and overbank areas. 


The discharges occurring each subdivision the cross section are de- 
termined follows: 


where 


total discharge cfs 


4 
\ 
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HYDRAULIC CHARACTERISTICS 
Section No. Main Channel 
Sub Total 
5702 1.27 1.173 3,020 3,020 
5707 263 4.99 2.920 38,000 38,000 
5711 7.21 3.732 83,000 83,000 
5713 547 8.17 110,000 110,000 
5713.5 570 8.55 118,000 118,678 
5714 8.77 125,000 129,150 
Overbank Area 
5713.5 0.25 0.397 678 
5714.0 192 0.62 0.727 4,150 
TABLE 


and 


Assuming 


dividing 


solving 


From 

definition 
curring 
The 
and the 
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FIGURE 


friction slope 


and corresponding elements for partial discharge and con- 
veyance capacity the subdivision under consideration. 


Assuming unit slope the above equations, they are written thus: 


Q=K,x1 and 


dividing one equation the other 


solving for 


From the preceding analysis, noted that the value, addition its 


definition conveyance capacity, can also designated the discharge oc- 
curring with slope equal unity. 


The friction head, determined averaging the slopes Sections 
and then multiplying the length illustrated below: 


4 
100 1000 100,000 1,000,000 
7 i 
q 
4 
j 


4 


1 2 

The velocity head, hy, derived weighting process using the partial 
discharges occurring each subdivision the cross section. This method 
identical the one used the Corps Engineers. (10) Velocities each 
segment are computed the basic equation, where the 


area the particular segment. The weighted velocity head results from the 
following relationship: 


Eddy losses are computed from the formula where 0.50 
for negative values and 0.10 for positive values. 
The following example application Method using the data 
fictitious stream, the Red Fox River, Colorado, for discharge 11,100 
cfs. The column numbers refer the columns Table 


Step 


Enter the cross-section designation under Column Column shows the 
assumed water surface elevation. For Section this has been established 
5714.0 from predeveloped rating curve. The main channel and overbank 
areas and conveyance values are entered Columns and 


Step 


The friction slope, computed from the equation 
shown Column Thus, the example, this equal 
0.00755. Note that the total (Column used. 


Step 


Column shows the subdivisional discharges computed the equation 


Substituting the values the main channel the example: 


similar computation made for the overbank area using its value 
4,150. 


Step 
The velocities Column are obtained dividing Column 


for each subdivision. Velocity the main channel equals 10,700 divided 
594. 


Step 


Column computed squaring Column and multiplying Column! 
The values each subdivision are totaled and entered. 


WATER SURFACE PROFILE COMPUTATIONS 


cfs 


Ascharge. = 20,000 


D 


Silver Fox River. 


SUBJECT: 
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Step 


The weighted velocity head, hy, computed from the following relation- 
ship: 


Column 
upstrean 


2gQ This 
Substituting the example: 
ers 
461 
which entered Column 12. the 
was done 
Step Again 
The distance between sections entered Column 
The 
elevation assumed and entered Column for Section Similar 
entries are made Columns and was done for Section described 
Step above. 
| 
Step 
The friction slope computed described under Step and entered 
Column Column shows the average the consecutive slopes follows: 
4 
0.00 
Column the friction head, equals Column multiplied Column 
Step 
Columns 10, 11, and are computed described Steps and 
however 
The algebraic difference velocity heads entered Column 13. the channel 
example, this equals 4.85 0.43 4.42. 
Step 
The eddy loss Column taken per cent Column for positive 
values per cent for negative values. This equals (0.10)(4.42) 0.44 mear 
the example. 
The 
Step the form 
Column the summation Columns and 14. Column equals 
Columns plus 13. 
Step Veloc 
Column obtained adding the water surface the downstream 
section (Column 17—Section the value shown Column 16. For 


this equals 5714.0 6.86 5720.86. When the result Column equal 


al 
ay, 


WATER SURFACE PROFILES 


Column the nearest tenth foot, the computations are begun the next 
upstream section. 


Method 


This method also trial and error procedure involving step compu- 
tations. used cases involving steady, nonuniform flow. However, 
differs from Method that the reach lengths representing the flow path be- 
tween sections are different for the main channel and overbank. The over- 
bank reach length could considerably shorter. Consideration also given 
the variability the hydraulic elements “n”, area, and hydraulic radius 
was done Method 
Again, the sketch Fig. representing the energy equation basic 
this method. Knowing the elevation Section the water surface elevation 
Section assumed and the energy equation checked. 
The determination the friction head, will discussed first. This 
accomplished applying Manning’s formula follows: 


q = an /3 g /2 


expressing where the length between sections. Substituting 


Substituting the equation 
h */e 
squaring and solving for 
2 
L ‘ 


however, since there are two more values for computed because 
the differences flow paths and conveyance capacities for the main 
channel and overbank areas, the summation used: 


mear friction head results from the equation where and 


are the friction heads computed for Sections and respectively. 


The discharges occurring each sectional subdivision are computed from 
the formula listed below: 


Velocity heads and eddy losses are treated the same Method 
application Method shown below using the data fictitious 


stream, Silver Fox River for discharge 50,000 cfs. Column numbers 
refer the columns Table 
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Step 
Enter the pertinent data the heading Table the example, these 


are the Silver Fox River water surface profile computations for discharge 


50,000 cfs. 
Step 

Enter the station number 209+40 under Column 
Step 


Under Column the main channel and overbank reach lengths between 
sections are entered. These are noted 2,260 and 1,955 feet between 
Sections 209+40 and 186+80 for the main channel and overbank reaches. 


Step 


The known water surface elevation 94.67 entered Column for 
Section 209+40. Areas each subdivision are entered Column 


Step 


values 973,000 and 604,000 are entered Column for the main 
channel and overbank area, respectively. 


Step 


Column determined dividing the individual values Column 


the square root Column For example, the main channel value comput 


follows: 


000 
20,500 
(2,260) 


Step 


The friction head, Column found squaring the quantity the 


total discharge divided the summation the values Column 
follows: 
2 


Step 


The subdivisional flows listed Column are computed multiplying 
the individual values Column the square root Column For ex- 


ample, the overbank discharge computed follows: 13,700 2.131/2 
20,000. 


Step 


Columns 10, and are determined shown Steps and under 


Method 


7 


Step 


The 
ginning 
surface 
sumed 
Column 
186+80 
entered 


Step 
Eddy 


Step 


Colu 
values 


Step 


Colu 
algebra 
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WATER SURFACE PROFILES 
Step 


The remainder the columns are omitted since Section 209+40 the be- 

ginning section. The next step, therefore, consists assuming the water 
surface elevation for the next upstream Section 186+80. This has been as- 

sumed 96.90 the example. Similar computations are made through 

Column described for Section 209+40. The velocity head, hg, for Section 

computes 0.45. The algebraic difference velocity heads 


Step 


Eddy losses Column are determined Step Method 


Step 


values Column shown the following computations: 


1.88 


Step 


Column equal the total Columns and 14. Column the 


Column the mean friction head loss, computed averaging the 

algebraic sum Columns and 15. 


the preceding water surface elevation (Section 209+40) Column 17. 
the example, this amounts 94.67 2.21 96.88. When balance the 
nearest tenth foot obtained between Columns and the computations 
are complete. 

The next step involves computation the friction head, between 
Sections 186+80 and 148+10 (next upstream section) using the lengths 3,870 
and 3,265 shown for the main channel and overbank. The water surface 
elevation Section 148+10 then assumed and new cycle computations 


begun using the same process outlined above. 


Step 
The water surface elevation Column calculated adding Column 
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TRANSITION FROM LAMINAR TURBULENT FLOW 


thorough review the paper which valid questions were posed concerning 
the interpretation the experimental results. 

The most important question the discussion concerned the propriety 
using the Blasius boundary layer the discussion stability the temporal- 


developing boundary layer pipe. Specifically, the question was raised 


whether not boundary-layer stability studies Tatsumi(2) would 
more applicable. Neither the Blasius boundary layer nor the Tatsumi bounda- 
layer the time developing boundary layer pipe, since the Blasius 
boundary layer applies steady-flow over flat plate without pressure 


gradient and the Tatsumi boundary layer applies steady flow pipe 


the inlet spite the differences origin, the Tatsumi boundary 


layer and the Blasius boundary layer are remarkably similar, shown 


Fig. which the distance from the boundary. Since the second deriva- 


tive with respect important stability analyses, velocity profiles 


may appear quite similar Fig. but result radically different 
neutral stability curves. Thus, further evidence the remarkable similarity 
between the Tatsumi and Blasius boundary layers shown Fig. which 
neutral stability curve for the Tatsumi boundary layer superposed upon 
the neutral stability curve the Blasius boundary layer. Therefore, there 
appears intrinsic merit using the Tatsumi boundary layer for 
qualitative discussion case for which the stability analysis not availa- 
ble. the other hand, the advantage using the Blasius boundary layer for 
comparison that the stability analyses are much more comprehensive than 
for any other boundary-layer velocity profile. 
The comparison the stability the temporally developing boundary laye 
pipe can also made with the stability the Polhausen boundary 
Fig. the plotted points were obtaineu from analog computer 
solution the Navier-Stokes equation with boundary conditions and physical 
constants Run 29. The analog computer solution for the velocity profiles 
shown Fig. the range boundary-layer thickness for which tran- 
sition turbulence was observed. The ratio the displacement thickness, 
the momentum thickness, about 2.43 for the four velocity profiles 
Shown Fig. This value corresponds Polhausen shape 
parameter Stability analysis this Polhausen boundary layer indicates 
that the critical minimum value the boundary-layer Reynolds number, 


Proc. Paper 1450, December, 1957, Carstens. 
Prof. Civ. Eng., Georgia Inst. Technology, Atlanta, Ga. 
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quired 
However, all the experimental runs, the inception The 
turbulence was always observed (see Table the paper 1450) value 
less than 3500. This observation tends substantiate the discusser’s 
statement “there some evidence that timewise developing flow may throug 
less stable than steady flow.” gradie 
spite the hazards using borrowed solution, that is, the stability sults 
the Blasius boundary layer, the general nature the boundary layer insta- 
bility the time developing boundary layer the experiments undoubtedly the 
identical that the borrowed solution. The fact that the borrowed stability 
analysis offered rational explanation the observed transition turbulence The 
indicative the author that the neutral curve stability the time de- The 
veloping boundary layer must very similar, both quantitatively and 
tatively, the neutral stability curve for the Blasius boundary layer. limite 
The discusser has correctly deduced the manner obtaining the “possible which 
neutral stability curves” which are shown Fig. 10, that is, the disturbance this 
wave length was chosen which would pass the “possible neutral stability until 
curve” through the experimental point having the highest transition Reynolds justm 
number. The hypothesis that the disturbance wave length, proportional 
the roughness height was based upon analogy with the Karman vortex appre 
trail behind circular cylinder. stable Karman vortex trail, the spacing burst 
between successive vortices proportional cylinder diameter and inde- bound 
pendent the flow variables. The analogy suggestive that the disturbance 
wave length, proportional only the roughness height or, this instal 
the step height the mismatched pipe boundaries the junction. The 


matching the pipe ends could not expected axisymmetric. Thus, 


| | 


known the author. any event the schematic drawing Fig. which 


many disturbance wave lengths, ranging downward magnitude from the 
maximum wave length the point greatest mismatching, could expected 
originate the pipe junction. Only the possible neutral stability curve for 
the longest disturbance wave length was plotted Fig. 10. The other shorter 
disturbance wave lengths originating the pipe junction would have neutral 


101 


stability curves lying below the plotted curve. The curve defining the limits 
neutral stability all disturbances would bounded vertical line 


Fig. 10. The above explanation was the basis for the statement “With this 
concept there conflict that some the points turbulence inception 


(indicated the circular symbols Fig. 10) are below the lower branch 
the curve shown.” 

The question was raised the discussion concerning the amplification 
disturbance which required produce turbulence inception. The stability 
analyses small disturbances are merely indications whether the dis- 
turbance will amplify decay. The stability analyses not indicate whether 
the amplification disturbance sufficient result transition turbu- 
lence. The presumption that, when disturbance reaches certain magni- 
tude, transition occurs. The order magaitude amplification required for 
transition not constant, implied the discusser, but must certainly 
depend upon the initial magnitude the disturbance. For example, infini- 
tesimal disturbance would require infinite amplification and finite dis- 
turbance, which equal greater than the limiting magnitude, would cause 
transition without amplification. Since the transition turbulerce was always 
observed some distance downstream from the pipe junction, the conclusion 


that the disturbances were small (but certainly finite) disturbances that re- 
quired some amplification before transition occurred. 


The discussion concerning Fig. appears result from misunderstand- 


ing. brief review the total transition occurrence will clarify the meaning 


this figure. Prior opening the downstream valve, the fluid rest 


throughout the pipe. the instant valve opening piezometric-pressure 


gradient applied the fluid the pipe. The unbalanced force system re- 


sults acceleration. The developing flow pattern can visualized 

central core with uniform velocity distribution and boundary layer which 
the velocity varies from zero the boundary the core velocity. The core 
velocity increases linearly with time during the initial stages the motion. 

The boundary layer zero thickness and increases with time. 


The flow laminar throughout. Meanwhile, small disturbances are being fed 
into the boundary layer downstream from the pipe junctions. There only 
limited range core velocity, and boundary-layer thickness, 6*, during 
which these disturbances will amplified shown Figs. and 10. During 
this period amplification, disturbance from the pipe junction amplified 
until turbulent burst formed. Once turbulence burst exists, large ad- 
justments the laminar velocity profile are effected. Since the pipe flow 
limited extent, these adjustments within the boundary layer also result 
appreciable velocity changes the central core adjacent the turbulent 
burst. Since the concept laminar central core riding over turbulent 
boundary layer appears untenable, the turbulence must penetrate into 
the central core some manner. Whether the velocity adjustment leads 
instability the central core whether the disturbance the turbulent 
burst sufficient cause instability and radial spread turbulence not 
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based upon study the film strip data and pressure data the downstream 
face the turbulent spot passed through the pipe outlet, appears 
rational picture the downstream face turbulent spot which moves 
toward the pipe outlet. Essentially, the same type downstream face 
turbulent spot reported Schubauer and Klebanoff(5) for flow over flat 
plate follows: 


“The turbulence is, fact, transported downstream with the free- 
stream velocity and the lag the surface due the time required 
for propagation inward through the slower moving air the laminar 
layer. Probably the chief significance the slower progress near 
the surface that gives rise overhanging leading edge.” 


Two errors the paper were pointed out the discussion. First, the 
number wave lengths between the point disturbance origin and the point 
turbulence inception would about for turbulent spot and about 
for turbulent spot based wave lengths 2.1 and 1.5 ro, re- 
spectively. Second, the X/D value for spot Run should rather 
than 112 Table question mark parenthesis should placed 
Table adjacent X/D value spot Run 21. 

The author thanks Professor Robertson for his carefully prepared dis- 
cussion the paper and thanks Messrs. Johnson and Meeks the Georgia 
Tech Analog Computer Section for their help obtaining the solutions shown 
Fig. 
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SOME EXPERIMENTS WITH EMERGENCY SIPHON 
Discussion Fred Blaisdell and Harold Humphreys 


FRED BLAISDELL,! ASCE and HAROLD HUMPHREYS,? 
addition the countries listed the author, has used large 
siphon spillways for the close control reservoir levels. But only India, 
the writers’ knowledge, has used both “saddle” and “volute” 
they call them. This discussion will confined largely the volute siphon 
and various modifications it. 

The volute siphon might defined morning glory spillway which has 
the entrance covered The dome greater diameter than the 
entrance and extends below the crest elevation the same manner does 
the top hood the saddle siphon. Also, like the saddle siphon, the volute 
siphon uses siphon breakers which run the top the dome. contrast 
the saddle siphon detailed Mr. McBirney, the volute siphon does not 
curve the lower leg upstream necessarily use device throw the nappe 
the crown, nor necessary submerge the exit insure priming. 

Indian literature credits Shri Ganesh lyar with the discovery the 
volute siphon 1933 and its experimental development Mandya 1935-36 
using siphon 2.5 feet diameter and feet Two volute siphons 
eight feet diameter with operating head feet were constructed 
Marconahally 1938 experimental measure. Eleven siphons feet 
diameter with drop from the design head-water level the invert the 
spillway exit feet were built 1947 Hirebhasgar. 

not necessary use the dome insure that the spillway will run full 
under suction. Examples this may found the Heart Butte, North 
Dakota, and Shade Hill, South Dakota, spillways the Bureau Recla- 
mation.(3) siphon defined spillway which requires that the liquid 
lifted above the headpool level, then the Dakota spillways not qualify 
siphons. They qualify siphons the definition is: 


closed conduit, portion which lies above the hydraulic grade line. 
This results pressure less than atmospheric that portion, and 
hence requires that vacuum created start flow. siphon util- 


izes atmospheric pressure effect increase the flow water 
through it. 


Proc. Paper 1807, October, 1958, Warren McBirney. 

Project Supervisor, Agri. Research Service, Dept. Agri., St. 
Anthony Falls Hydr. Lab., Minneapolis, Minn. 

Head, Hydr. Research Section, State Water Survey, Urbana, 
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The Dakota spillways qualify siphons because, under full flow, the vertical 
legs are under pressures less than atmospheric and lie above the hydraulic 
grade line. 

The possibility that the pressures within these spillways may low 
cause cavitation should recognized. The Hirebhasgar volute siphons(1) 
India are examples. These siphons consist dome-covered morning 
glory entrance feet inches diameter which converges drop inlet 
shaft feet diameter, degree elbow, and horizontal leg feet 
diameter. The drop level between the crest and the exit center line about 
feet—considerably excess absolute zero pressure head even consider- 
ing hydraulic losses. Cavitation damage should have been expected and did 
Cavitation damage also has occurred the Marconahally siphons. 
This should have been anticipated since one-sixteenth size model tests indi- 
cate prototype pressures below absolute zero over much the drop inlet and 
inside the elbow.(5) 

the Dakota spillways the Bureau Reclamation, the possibility 
extremely low pressures was recognized and means were taken prevent 
their occurrence. the case the Heart Butte spillway, lip used 
cause the stream break away from the wall the drop inlet. This insures 
atmospheric pressure that point and that negative pressures the drop 
inlet will not exceed predetermined magnitude. The drop through the Shade 
Hill spillway does not appear sufficient cause cavitation problems. 

Full flow can achieved for spillways the morning glory type without 
special priming devices the inlet properly formed. The writers are 
the opinion that the Hungry Horse spillway(3) would flow full the head over 
the crest raised high enough. this were occur, serious cavitation 
would likelihood, view the 487 foot drop through the spillway, unless 
means were taken prevent it. Presumably the maximum head over the 
crest sufficiently low and the barrel sufficiently large insure part full 
flow all times. important that designers recognize that morning glory 
spillways can flow full, order know what design against they want 
part full flow, know what design for they want full flow. 

The Dakota spillways required priming depriming aids other than 
shaping the inlet properly. The head-discharge curves for the Heart Butte 
spillway presented Peterka(6) show that the flow through the spillway 
controlled weir low heads. high heads, the control pipe flow 
when the inlet submerged and the drop inlet full. the higher weir 
heads, the falling.water must partially seal off the drop inlet and the increased 
head must create demand for more water than can supplied over the 
crest weir that headpool level. The excess demand supplied air en- 
trainment and air flow through the spillway. This is, the writers believe, the 
partialization for which Mr. McBirney looking. Partialization achieved 
here without special devices. 

Only relatively large spillways have been discussed this point. Large 
closed conduit spillways designed flow full, like the Heart Butte and Shade 
Hill structures, are few number. However, the Soil Conservation 
Service (SCS) began building this type spillway relatively smaller sizes 
the 1930’s. Closed conduit self-priming spillways built annually under 
Public Law No. 46, 74th Congress, entitled “An Act Provide for the Pro- 
tection Land Resources Against Soil Erosion, and for other Purposes;” and 
Public Law 566, 83rd Congress, known the Watershed Protection and Flood 
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Prevention Act, number the and many them operate de- 
sign capacity each year. The hydraulic design these SCS spillways 
identical the design methods which should used for the relatively larger 
spillways the morning glory type. The senior writer has outlined the theory 
previously(8) and recently has authored series describing both 
laboratory and field tests number different inlet configurations. 
may interest mention passing that the hydraulics involved that 
taught elementary courses but that the application requires thorough 
knowledge the conditions under which each the elementary hydraulic laws 
are applied. 
great many the SCS-designed closed conduit spillways consist 
vertical drop inlet from the base which barrel carries the flow through 
the dam. The slope the barrel may either steeper flatter than that 
the hydraulic grade line. the drop inlet properly proportioned, the spill- 
way will flow full whatever the barrel slope may even the outlet not 
submerged. Pressures within the barrel will less than atmospheric and 
the barrel will lie above the hydraulic grade line the barrel slope steeper 
than the hydraulic grade line. The spillway then acts siphon the broad 
definition presented earlier accepted. 
series tests drop inlet closed conduit spillways that are pertinent 
this discussion was initiated 1957 the Urbana Hydraulic Laboratory 
the State Water Survey under the direction the junior writer. This 
spillway consists circular drop inlet five barrel diameters (5D) high 
5/3 barrel diameters (5D/3) diameter. The barrel diameter three 
inches. The per cent slope which the barrel has been placed initially 
considerably steeper than the hydraulic grade line slope. Water enters the 
drop inlet from the entire periphery the crest. The head-discharge curve 
for this spillway has the form shown Fig. Section the head-discharge 
curve represents weir control with the weir the crest the drop inlet. 
Section also weir control but air carried through the spillway with the 
water. Section represents pipe control. The head-discharge curve and the 
flow have been fully described elsewhere. (10) 
The primary purpose the State Water Survey study de- 
termine the performance flat plate anti-vortex device supported above 
the crest and develop design rules for the form and dimensions the plate. 
The initial flat plate two drop inlet (riser) diameters Dp) diameter 
with three guide vanes, 120° apart, attached the plate. The guide vanes ex- 
tend radially outward from the outside edge the drop inlet. The solid curve 
Fig. for the case where the anti-vortex plate located 
higher above the inlet crest, that is, within the range pipe control repre- 
sented Section When the plate was lowered into the weir control ranges 
was found that the plate sealed off the inlet soon the water 
level reached the bottom edge the plate. This prevented access air 
the drop inlet and the flow caused the creation partial vacuum. The 
vacuum sucked air under the plate sufficient quantity satisfy the demand. 
Increasing the water flow caused little change headpool level until 
the spillway flowed full water alone. Air was sucked into the spillway 
amounts maintain vacuum just sufficient permit the water flow take 
place constant headpool levels. The dash lines Fig. show the head- 
discharge curves for successively lower elevations the anti-vortex plate. 
The headpool level essentially steady, exhibiting none the fluctuations 
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shown Mr. McBirney’s Fig. this not the partialization and steady 
headpool level that Mr. McBirney looking for? The possibility positively 
preselecting the priming head within certain presently undefined limits 
setting the flat plate the correct height above the crest should not over- 

give actual figures for comparison, Mr. McBirney’s spillway apparently 
has throat area 0.5 0.5 square feet. The Dlinois spillway has 
barrel area (0.25)2 0.0491 square feet. The linear multiplication 
factor secure comparable priming heads then 3.19. Plate 
heights low 5/8 inch caused the spillway prime. Increasing 
this the scale Mr. McBirney’s tests, the priming head becomes 0.17 feet. 
This corresponds Mr. McBirney’s full flow priming head 0.52 feet for 
the standard siphon design and 0.42 feet for the proposed design. 

Proper positioning the flat plate anti-vortex plate used siphon 
primer results single-valued head-discharge curve. The curves two 
loops shown Mr. McBirney Figs. and for the standard and proposed 
designs not exist for the flat plate priming device the positions shown 
Fig. However, two loops can obtained for the drop inlet entrance ex- 
tending skirts from the periphery the cover plate down below the crest 
the weir. The spillway would then prime rising stages the same head 
the skirts had not been used, but the prime would held until the pool level 
reached the bottom the skirts, the siphonbreaker one used. The 
skirted arrangement is, course, what called vortex siphon spillway 
India. 

The priming time, used Mr. McBirney for the saddle siphons, does 
not exist for the drop inlet closed conduit spillway. Priming for the latter 
spillway takes place gradually the headpool level rises and does not occur 
relatively suddenly for the saddle siphon. 

The gradual increase and decrease flow obtained with the drop inlet 
siphon also reduces the erosion damage downstream from the outlet which 
Mr. McBirney says results many locations from the prime-break-prime 
cycle intermittent operation. 

The writers not say that the drop inlet closed conduit spillway, either 
without with the flat plate anti-vortex device and primer, complete 
solution all siphon problems. They maintain, however, that the potenti- 
alities the drop inlet siphon should not overlooked. 
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SYNTHETIC FLOOD FREQUENCY4 


Discussion Max Kohler 


synthetic, flood frequencies can best described elaboration and im- 
provement the so-called method.” involves the assumption that 
the n-year flood can derived from the n-year rainfall even though the two 
ory events not necessarily bear the relation cause and effect. 
Specifically, the proposed method involves (1) determination the n-year 
rainfall for the basin with specified area and concentration time (2) con- 
version n-year rainfall n-year runoff volume and (3) determination 
the flood peak corresponding the runoff volume. This discussion concerns 
primarily steps (1) and (2) the procedure since believed that Snyder’s 
any several other available techniques are adequate for converting run- 
off volume corresponding peak flow. 
Before taking specific aspects Snyder’s procedure, perhaps best 
consider the basic underlying assumption that simple relation exists be- 
tween the n-year flood peak and the n-year rainfall. From probability point 
view, this implies that the composite probability for all factors which af- 
fect the relation between rainfall amount and peak flow unity. Therefore, 
the only direct means relation n-year rainfall and flood peaks through 
statistical analysis n-year data. This means that flood-frequency data 
must used deriving the generalized relations. There indication 
the part the author that flood-frequency data are necessary were used. 
Certainly there was some basis for deciding that the required relation should 
give “slightly higher runoff than the average” (indicated average monthly 
rainfall and runoff). Wouldn’t one expect different rainfall-runoff relation 
for small basins hrs) where the major floods result from summer 
thunderstorms than for the larger basins where such floods result from gener- 
rains wet soil? believed the author has tended account for this 
variation unrealistic reduction point rainfall for and drainage 
area. Thus chose derive area-reduction relation from storm data 
using one properly derived from frequency data.(1) not 
clear why his results indicate flatter curve than one based frequency 
data, while just the opposite would expected. 

visualized that Snyder’s procedure constitutes method generaliz- 
ing, regionalizing, the basis existing flood-frequency data, certain 
modifications would appear appropriate. Although not absolutely neces- 
suggested that the empirical tie between n-year rainfall and flood- 
peak confined the rainfall-runoff relation. this case, the n-year 


Proc. Paper 1808, October, 1958, Franklin Snyder. 
Chf. Research Hydrologist, Weather Bureau, Washington, 
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runoff can computed from Snyder’s relations for all basins with established first 
flood-frequency curves. These values can then related corresponding 
values n-year basin rainfall derived from the point-rainfall frequency re- 
lation (author’s Fig. and one areal rainfall duration and point and 
rainfall.(1) This the manner which the accompanying chart was de- 
veloped. exac 
The values peak discharge (Qp) for various basins and return periods the 
were obtained Gumbel analysis with appropriate adjustment give partial- was 
duration data. The average rate direct runoff for time was obtained two 
use the author’s formula, Qp/ 500A, where the basin area labe 
square miles. 
The first part the accompanying relation, namely, the area curves the 
upper left, was derived from the 10-year rainfall curve the author’s Fig. 
The 10-year rainfall amounts for various concentration times, durations, 
were first plotted and the zero-area curve was drawn through the points. The 
curves for other sizes area were then constructed the basis the depth- 
area relation(1) given Weather Bureau Technical Paper 29. Although tens 
the depth-area relation based rather limited data, considered 
the best available for the purpose. The second part the writer’s relation, yeal 
the return-period curves the lower left, was obtained entering the rela 
: TEN- YEAR RAINFALL INTENSITY (1) IN INCHES PER HOUR and 
Tp (yr) 
NW Br Anacostia R. or Colesville, Md 21.3 3 138 261 407 
a Appomattox R. at Farmville, Va 306 31 47 84 
Linganore Ck. nr Frederick, Md 82.3 7 65 96 142 
: 4 Monocacy R at Jug Bridge, Md 817 24 55 860 «131 3. 
< Octararo Ck. nr Rising Sun, Md 193 9 72 151 264 
8 Rappahannock R. nr Fredericksburg, me. 1599 24 48 94 160 7 
5 2 Seneca Ck at Dawsonville. Md 101 7 4 109 189 
é Slate R or Arvonia, Va 235 31 42 72 116 


T)-YR RAINFALL INTENSITY 


0.02 0.04 0.06 oO. 02 O06 | 2 


ToYR RUNOFF INTENSITY (I,) IN INCHES PER HOUR 
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first quadrant with going horizontally the zero-area curve, then down- 
ward the appropriate rainfall intensity (indicated the author’s Fig. 2). 
Each point determined was labelled with the corresponding return period 
and smooth family curves were fitted the points. (2) 

should pointed out that the return-period curves cannot drawn 
exactly fit areas all sizes, but the errors are insignificant with respect 
the errors the procedure whole. The final curve family curves) 
was derived plotting n-year basin rainfall intensity computed from the first 
two quadrants against corresponding values n-year runoff intensity and 
labeling the points with 

The accompanying table presents the basic data for eleven natural basins 
the general vicinity Washington, C., which were used deriving the 
These basins were selected because flood-frequency curves based 
long records observed and synthesized(3) annual peaks were available 
were unit hydrographs for estimating 

The application the writer’s relation for estimating n-year runoff in- 
tensity straightforward and needs further amplification here. The pro- 
cedures described Snyder for obtaining and for converting the 
year peak flow are equally applicable when using the writer’s accompanying 
relation. 

conclusion, the writer wishes commend the author for presenting new 
and stimulating ideas important phase hydrologic design. 


“Rainfall Regime, Part The Ohio Valley,” 
Weather Bureau Technical Paper No. 29, June 1957. 


Paulhus, pp. 311-321, McGraw-Hill, 1958. 


“Flood Frequencies Derived from Rainfall Data,” Paulhus and 


Miller, Proc. Paper 1451, Jour. Hydraulics Div., ASCE, Vol. 83, 
December 1957. 
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DIVINING RODS VERSUS HYDROLOGIC DATA AND 


Discussions Erhard Dittbrenner and Yevdjevich 


ERHARD DITTBRENNER, ASCE.—Mr. Langbein presents scholar- 
and comprehensive review the status the most critical need hy- 
drology. would difficult add anything useful his presentation hy- 
drology. However, when includes “waterwitching” hydrology—even 
only strawman knocked down and used reference—it time 
separate fact from fancy. 

Mr. Langbein states: “Modern hydrology offers support water 
witchery.” certain limited sense this statement may accurate enough. 
His further question: can scientific about witchery,” obviously 
rhetorical, his implication that since witchery perhaps the realm 
fakery, cannot scientific. 

The answer his question is, however, unqualified “Yes.” can 
scientific about “witchery.” so, however, must that portion 
physical science which has with electrical phenomena. This might seem 
take out that portion science which have artifically labeled “hy- 
drology.” not important argue the point, but seems pertinent 
point out that push back the frontiers our knowledge our artificially 
prescribed categories science becomes clear that Nature has such 
arbitrary divisions. More and more evidence substantiates belief perhaps 
more ancient than water witching—that there oneness unity the 
Universe. 

The writer happens “water witch” dowzer. Many years ago 
was invited try his hand young dowzer who suspected skepticism the 
art attitude. The dowzer pointed out that many people have the dowzer’s 
capabilities without realizing it. One trial demonstrated the truth that 
statement and the basic reason why water witchery has persisted over the 
world and through the ages. Once person experiences the phenomena which 
turns the switch down, will never bothered the scoffing the “scien- 
mind. the writer’s case, the switch broke his hands his first 
trial. has never done since, but since never worked the art very 
much this may not very significant. 

Since there had some scientific explanation for the phenomena, search 
literature and query among those more qualified science than himself 
finally brought some answers. The “dry” holes can explained the 
knowledge that the forces involved the phenomena are obviously too small 

break the twig, but act certain muscles which are thereby impelled 
bring about the action the switch without volition the part the dowzer. 


Proc. Paper 1809, October, 1958, Langbein. 
Hydr. Engr., Bur. Public Rds., Albany, 
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The forces involved vary constantly—during the day, with the weather, and 
from season season. Some experienced dowzers are aware this. Others 
obviously are not. Since the dowzers are not aware the source the phe- 
nomena their powers, they are unable explain these variations. 

Two scientists England have studied this phenomena extensively. Maby 
and Franklin, “The Physics the Divining Rod,” London, 1939, cover the 
subject comprehensively, but means exhaustively. The phenomena has 
been measured. Instruments have been devised determine quantitatively 
the field which the dowzer works make his determinations. People have 
been checked determine their powers “divination.” About one person 
ten has these powers some degree. Their investigation did not extend 
the point determining the reasons why some and others not have 
them. far this writer knows, one has made any effort this di- 
rection—in fact, since Maby and Franklin’s publication, further scientific 
literature that nature the subject has appeared the writer’s knowledge. 

the writer’s own experiences, routine laboratory experiment during 
college required the measurement each student’s bodily electrical re- 
sistance. His resistance was way out line with that any other person 
the group. The notes are gone but the memory says was lower—perhaps 
third fourth the “normal.” Some interim reading—also not then noted 
that persons who become alcoholics must more careful around 
electricity than others, they are more easily electrocuted. There 
implication that persons with low electrical resistance tend toward become 
alcoholic victims, but there implication that bodily resistance electric 
currents vary greatly and these variations have significant physical and 
physiological results. Medical research some years ago, has es- 
tablished the fact that the human body generator well conductor 
electricity. logical that generation and restivity would bring about 
difference the electrical field about the human body and its reaction with 
the field about water “veins.” 

Here, physiology and the physical sciences appear have common inter- 
est. far, nothing has gone beyond what Maby and Franklin have reported, 
far the writer knows. would seem interesting field for both 
the medical and the earth scientist explore. 

The question now arises—“so what?” Granted that the dowzer has definite 
powers, useful the geologist the hydrologist? Mr. Langbein 
himself answers that question. the absence completely detailed geologi- 
cal and hydrologic subsurface knowledge, the dowzer certainly can useful. 
the writer’s opinion, the likelihood any other knowledge supplanting him 
appears remote. far the writer’s knowledge geology and hydrology 
permit, neither the latter branches science will produce the detailed 
knowledge available the informed and experienced dowzer. seems possi- 
ble that extension the Maby and Franklin research could enhance the 
dowzer’s powers, and mayhap even elevate him “professional” category. 

This would appear subject useful and profitable research. There 
certainly further basis for scoffing. 


7 


ASCE.—Two very interesting problems are dis- 
cussed Mr. Langbein: contrast the divining rods 


Hydr. Engr., Foreign senior scientist, National Bureau Standards, 
Washington, 
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hydrologic decision-maker versus the modern basic-data programs; con- 
thers trast the increase gathering basic-data versus the deficiencies re- 
search necessary guide, analyze and interpret the basic data, and obtain 
both the general principles and the regional hydrologic characteristics. 
laby clear that modern hydrology does not support “forked stick” pro- 
the cedure which does not have serious scientific background. But the hydrolo- 
gist must aware some facts which contribute the persistence the 
dowzer’s activities, that the hydrologist may overcome his seemingly ap- 
parent disadvantages. Ground water hydrology, frequently very complex, im- 
poses some uncertainty decisions upon the geohydrologic specialist, while 
the dowser’s short cut way decisions makes spectacular impression. The 
dowser often man well acquainted with local conditions, and makes de- 
cisions analogy, based experiences. sometimes kind living 
ific map, easy for farmer consult. available, and what important, 
cheap. These are the reasons why the use the divining rod persistent 
and will continue for relatively smaller activities until replaced 
good engineering procedures, competitive from other points view. 
The collection basic data hydrology has been started and made under 
some general assumptions for the expected water resources development and 
control. The basic concepts laid for the collection hydrologic data depend 
upon the future water resources development. The first concepts were based 
the expected individual and mostly single purpose projects. Next step 
development was the multiple purpose use water resources. The future 
will probably ask the total and very integrated water resources development 
state great demands for water control and use. The basic concepts 
governing the collecting hydrologic data today should anticipate the situation 
tor water resources development decades more from now. This 
another reason support Mr. Langbein’s emphasis for better scientific 
ith foundation general collection and interpretation basic-data. 

When demand for higher scientific efforts hydrology forwarded, the 
inter- problems connected with have emphasized. The efficiency 
ted the research hydrology important question. Surveying the hydro- 
logic literature the countries with relatively abundant funds and large 

number people employed research, the impression that the results 
finite not correspond the efforts. The quality the results and the possibilities 

for the useful applications are rather limited comparison with the quantity 
published papers and other literature. The greatest results for given funds 
and their best application the current needs and problems can achieved 
only timely and thorough preparations for expending research. 
ogy The timely solution many problems can help the rapid development 
research hydrology, for example: the systematic preparation highly 
educated engineers and scientists hydrology; the development good pro- 
cedures putting forward and selecting the topics studied; the ex- 
change experiences, ideas and new methods; the type work, teams, 


individual both; etc. 
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QUEUING THEORY AND WATER 


GOULD.!—The author’s paper shows that the techniques mathe- 
matical statistics may applied the solution the problem determining 
required water storage, and also estimating the probable behaviour 
reservoir already constructed. 

The service function: 


D=b+kS 


may interpreted several different ways, each which gives different 
result. 


avoid confusion, the writer has adopted the following notation: 


Dx, Average rate discharge the interval between time and time 


Ix, Average rate inflow the interval between time and time 


The unit storage used the writer the same that used the 
author, being equal unit flow for unit time. 

The first interpretation the service function considered that the 
discharge during interval time determined the storage the end 
the interval. This type control would difficult achieve practice, 
but considered because gives the same solution that obtained the 


author. 
and similarly 


The basic equations obtained the author may verified subtraction, and 
substitution 


Another interpretation the service function could that the instantane- 
ous discharge dependent the instantaneous value storage, i.e. 


Proc. Paper 1811, October, 1958, Langbein. 
Engr., Horsham Waterworks Trust, Vic., Aust., formerly lecturer Civ. 
Eng., Univ. New South Wales, Aust. 
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subtraction and substitution, this gives 


the time interval sufficiently small, Dx, x,1 approximately equal 
1/2 Substituting this value and rearranging gives: 


instead the author’s 


yet another case, may taken that the discharge could determined 
the amount water storage the beginning interval. This type 
control would suitable for irrigation area containing both semi- 
permanent (e.g. citrus vines) and annual (e.g. market garden) development, 
where certain minimum amount water required maintain the semi- 
permanent vegetation, whilst the area annuals planted any year may 


varied according the amount stored water available the beginning 
the year. This gives the equations: 


and 


Subtracting, and substituting for Sx-1) gives 


x=1 9x 


estimate the discharge the nth year, given successive substi- 
tution “D” values, is: 


From this, the standard deviation the estimate the discharge can 
calculated as. 


instead 
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DISCUSSION 


given the author Eq. (6). 
comparison between the author’s Eqs. (12), (13), and (14), and the corre- 
sponding equations derived the writer for the case where the discharge 
during period based the storage available the beginning the peri- 
od, given the following table: 
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Ean. Discharge based Discharge based 
No. final storage. initial storage. 


nent, 
mi- 

may 
ing 


noted that where the discharge determined the amount 
stored water the beginning the period under consideration, the equations 
give when this case, any flow excess the minimum flow 
any year may considered being stored for use the immediately 
following year, thus giving time make plans for its utilization. Under such 
conditions, the storage required provide 100 per cent utilization would 
the difference between the maximum probable inflow, and the minimum proba- 
ble natural inflow. This verified the equation which gives the required 
storage 

the case considered the author, the storage required increases 
rapidly approaches and becomes infinite when 

The accuracy the “probability routing” methods described the author 


examples and the paper, limited only the number storage 
intervals used, the number successive approximations made, and the ever- 
present errors sampling. Electronic computation using standard 
programmes for the solution linear simultaneous equations has been used 
+ 


the writer, who found that smaller intervals storage volume could 


used, and that the need for successive approximations was eliminated because 
the solution was provided direct manner. 
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DIVISION ACTIVITIES 
HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


NEWS 


April, 1959 


PURPOSE THE HYDRAULICS DIVISION 
(Quoted from the Official Register) 


“The advancement and dissemination knowledge relating the oc- 
currence water nature and its behavior structures, water courses, 
and underground. 

“In particular the field the Hydraulics Division shall embrace meteor- 
ology and hydrology they affect the engineer, fluid mechanics engineering 
usage, and applied hydraulics branch engineering science which 
furnishes the basis for hydraulic design and for the practical use water 
the different specialized branches hydraulic engineering.” 


COMMITTEES SPONSOR FOUR SESSIONS CLEVELAND 


Four committees the Division have each prepared session for the 
ASCE Convention Cleveland May 4-8, 1959. The four sessions were 
published tentative the February Newsletter, and the latest information 
indicates that changes have been made the programs. 


HYDRAULICS CONFERENCE FORT COLLINS, COLORADO 
July 1-3, 1959 


The eighth national Conference the Hydraulics Division will held 
Fort Collins, Colorado July The hosts this year are the Colorado 
and Wyoming Sections ASCE and Colorado State University. There will 
three days technical sessions and one day extensive field trip through 
the beautiful Colorado Rockies. addition the professional sessions there 
will social sessions which fun and good fellowship can enjoyed all. 


Conference Headquarters Colorado State University Campus 


Aylesworth Hall, beautiful new dormitory the University Campus, will 
the housing center for the conference. This close the new Engineering 


Note: No. 1959-16 part the copyrighted Journal the Hydraulics Division, 
Proceedings the American Society Civil Vol. 85, April, 1959. 


Copyright 1959 the American Society Civil Engineers. 
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Auditorium where all conference sessions will held. Rooms will availa- 
ble the following special rates: $2.50 single, $2.00 per person double. Ex- 
cept for the special events, meals will available the following prices: 


Adults Children 
Breakfast $0.75 $0.65 
Lunch 1.00 0.95 
Dinner 1.25 1.10 


addition the facilities listed above, there are several hotels and 


motels available. Write the Housing Committee alternate accommodations 
are preferred. 


Facilities for Technical Sessions and Tours 


The auditorium the new Engineering Center has been reserved for the 
technical sessions. conducted tour the Hydraulics Laboratory and other 
engineering facilities one the features the program. 

guided tour the Big Thompson Irrigation and Power Project has been 
planned. This trip includes ride across the Continental Divide over Trail 
Ridge Road (elevation 12,000 ft) scenic Rocky Mountain National Park. 
trip through the Hydraulics Laboratory the Bureau Reclamation 
Denver also the schedule. 


Social Events 


After the field trip Tuesday, June 30, social hour with cocktails has 
been planned the evening. barbecue will held Wednesday evening 
and special banquet Thursday evening. All these events are open for 
family participation. 


Additional activities are being planned for the ladies and children at- 
tendance the conference. 


Vacation Possibilities 


Fort Collins located near the heart scenic Colorado. Bring the family 
and plan spend vacation the Rockies where the scenery magnificent, 
the fishing good, and many the housing facilities include swimming pools 
and other recreational facilities. Hiking, horseback riding, water skiing, 
tennis, and other facilities are available for families who enjoy contact with 
the great outdoors. Since normal tourist are heavy during the 
summer months, arrangements for vacation the mountains should made 
early through the Housing Committee the Conference. 


Six Committees Sponsor Technical Sessions 


The tentative program for the Hydraulics Division conference follows: 


Sponsored the Hydraulic Structures Committee 
Presiding: Kindsvater, Chairman, Executive Committee, and 

Peterka, Member, Hydraulic Structures Committee 
THE MANIFOLD STILLING BASIN, Maurice Albertson, Director, Re- 
search Foundation, Colorado State University, Fort Collins, Colo., and Gene 


Fiala, Highway Design Engineer, Bureau Public Roads, Portland, 
Ore. 
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THE VORTEX CHAMBER AUTOMATIC FLOW-CONTROL DEVICE, 
Richard Kolf, Assistant Professor, and Paul Zielinski, Assistant In- 
structor, Marquette University, Milwaukee, Wis. 

NEW CONCEPTS TUNNEL SPILLWAY DEFLECTOR-BUCKET DESIGN, 
Rhone and Peterka, Engineer, Bureau Reclamation, Denver, 
Colo. 


Sponsored the Hydromechanics Committee 
Presiding: Arthur Ippen, Vice-Chairman, Executive Committee, and 
Donald Harleman, Member, Hydromechanics Committee 


SYMPOSIUM HYDRAULIC ONE 


UNIFIED CONCEPT DYNAMIC SIMILARITY FLUID MODELS, 
Donald Harleman, Associate Professor Hydraulics, Massachusetts 
Institute Technology, Cambridge, Mass. 

SIGNIFICANCE AND APPLICATION FROUDE AND REYNOLDS 
NUMBERS CRITERIA FOR SIMILITUDE, Liu, Associate Civil Engi- 
neer, and Albertson, Director, Research Foundation, Colorado State 
University, Fort Collins, Colo. 

DYNAMIC SIMILARITY WAVE AND TIDAL MODELS, Gerritsen, 
Associate Professor, Coastal Engineering Laboratory, University Florida, 
Gainesville, Fla. 

SEPARATION AND SCALING INERTIA AND VISCOUS FORCES 
PERIODIC FLOWS, Gershon Kulin, Hydraulic Engineer, Fluid Mechanics 
Section, National Bureau Standards, Washington, 


Sponsored the Tidal Hydraulics Committee 
Presiding: Tiffany, Secretary, Executive Committee, and 
Disney, Chairman, Tidal Hydraulics Committee 


SYMPOSIUM INSTRUMENTATION 


NEW INSTRUMENTS DEVELOPED THE U.S.G.S. FOR THE MEASURE- 
MENT TIDAL FLOW, Barron, Hydraulic Engineer, Research 
Section, Surface Water Branch, Geological Survey, Columbus, Ohio. 

RECENT DEVELOPMENTS THE COAST AND GEODETIC SURVEY 
EQUIPMENT FOR THE MEASUREMENT AND REMOTE RECORDING 
TIDES AND CURRENTS, Goodheart, Division Tides and Currents, 

Coast and Geodetic Survey. 

PRESSURE-TYPE TIDE AND WAVE RECORDER, Bowler, Instru- 
ments and Control Systems Laboratory, National Research Council Canada. 
SPECIAL DEVELOPMENTS OCEANOGRAPHIC INSTRUMENTS, 
Snodgrass, Head, Special Developments Division, Scripps Institution 

Oceanography. 


Sponsored the Sedimentation Committee 
Presiding: Harold Martin, Member, Executive Committee, and 
Carlson, Member, Sedimentation Committee 


SEDIMENTATION AND CONSERVATION WESTERN ARID LANDS, 


Peterson, Branch Area Chief, General Hydrology Branch, Division Water 
Resources, Geological Survey. 

THE USE ULTRASONICS THE MEASUREMENT SUSPENDED 
SEDIMENT SIZE DISTRIBUTION AND CONCENTRATION, Gordon 
Flammer, Assistant Professor Civil and Irrigation Engineering, Utah 
State University. 
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Sponsored the Hydrology Committee 
Presiding: Maurice Dickinson, Member, Executive Committee, and 
Koelzer, Member, Hydrology Committee 


SYMPOSIUM THE USE ELECTRONIC COMPUTERS 
HYDROLOGY 


APPLICATION ELECTRONIC COMPUTERS THE SOLUTION 
HYDROLOGIC PROBLEMS THE BUREAU RECLAMATION, Francis 
Swain and Riesbol, Bureau Reclamation, Denver, Colo. 

THE ROLE ELECTRONIC COMPUTERS HYDROLOGIC STUDIES 
THE TVA, Snyder, Tennessee Valley Authority, Knoxville, Tenn. 

COLUMBUS RIVER SYSTEM POWER ANALYSIS USE DIGITAL 
COMPUTER, David Lewis and Loren Shoemaker, Corps Engineers, 
Army, North Pacific Division, Portland, Ore. 

THE USE HIGH-SPEED COMPUTING EQUIPMENT WATER RE- 
SOURCES INVESTIGATIONS THE GEOLOGICAL SURVEY, 

Harbeck, Jr., and Isherwood, Geological Survey, Denver, Colo. 


Sponsored the Flood Control Committee 
Presiding: Chamberlain, Co-Chairman, Fort Collins Hydraulics 
Conference, and Arno Lenz, Member, Flood Control Committee 

FLOOD CONTROL PROBLEMS THE SOUTH PLATTE RIVER BASIN, 
COLORADO, Stanley Miller, Supervisory Water Conservation Project Engi- 
neer, Colorado Water Conservation Board, Denver, Colo. 

CHANGES URBAN OCCUPANCE FLOOD PLAINS THE UNITED 
STATES, Gilbert White, Chairman, Department Geography, University 
Chicago, Chicago, 

NEW APPROACH LOCAL FLOOD PROBLEMS, Gen. Herbert 
Vogel, Chairman the Board, Tennessee Valley Authority, Knoxville, Tenn. 

SUGGESTED OUTLINE LEGISLATION ACCOMPLISH FLOOD 
PLAIN ZONING, Author announced. 


| 


HYDRAULICS DIVISION MEMBERS SERVE THE 
INTERNATIONAL STANDARDS ORGANIZATION 


Announcement has been made ASCE Headquarters the appointments 
Messrs. Rolland Carter, Frank Campbell, and Norbert Leupold 
ASCE representatives Subcommittee ISO/TC 30/SC1 the International 
Standards Organization, with subject: Flow Measurements Open Channels. 
Mr. Carter will serve with Work Group No. which deals with velocity- 
area method measurements; Mr. Leupold with Work Group No. will 
concerned with notches, weirs and flumes; and Mr. Campbell will serve with 
Work Group No. regarding glossary terms. 
meeting the Subcommittee was held London, 8-15, 1958, 
attended Mr. Carter who was the representative the Hydraulics Division, 
ASCE and the American Standards Association. 


REVIEWS TECHNICAL PAPERS ARE INADEQUATE 


The Hydraulics Division during the past three four years has seen 
gradual increase the number papers received for review and consider- 
ation for publication. This increase has reached the point where adequate 
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reviews members the Division are not accomplished great many 
instances. the result number papers are recommended for publi- 
cation without critical review determine their real importance. The 
problem stems not only from the increase the number papers but also 
from lack information who most qualified review paper and 
the time required for the reviewing process. Also, the problem further 
complicated when the reviewer and/or the author are not fully aware the 
standards and procedures established for processing technical papers. 

The Division Executive Committee its meeting last November appointed 
task force composed the present and past chairman the Division’s 
Publication Committee (Messrs. Ball, Fortson, Jr., and 
Dewey, Jr.) for the specific purpose studying the technical paper review 
procedures with view improving and accelerating the process. The suc- 
cess their efforts must ultimately depend upon understanding the 
publication procedures both authors and reviewers and increased inter- 
est resolving the problem all members the Division. can start 
“We Can Improve the Processing Technical Papers,” and com- 
plete the form the inside back cover this Newsletter. 


CAN IMPROVE THE PROCESSING TECHNICAL PAPERS 


you know how much work goes into the processing your technical 
papers before they appear the Journal? doubt you have pondered this 
question and wondered about the answer. The following brief outline and 
information are offered give you least partial answer the question, 
and provide you the opportunity help speeding and improving the 
process: 


Two copies your paper are sent from the New York Office the 
Publications Chairman the Division for review, that is, the re- 
quired three copies are sent the New York Office, instructed 
the inside front cover each issue the Journal. 

The two copies received the Publications Chairman are sent two 
reviewers, experts the field covered the paper, insofar the 
Chairman knows. The reviewers are selected from list 
names collected over the years the chairmen the Committee. Ad- 
ditional competent reviewers are needed. 

The reviewers make detailed study the paper its technical 
contents, adequacy and suitability for publication the Journal and re- 
port their findings the Chairman. 

From the reviewers’ reports the Chairman decides whether the paper 
reviewed third time, published, declined. least two 
unfavorable reports are required decline publication. The paper 
either accepted presented, accepted condition that changes are 
made clarify satisfy questions and requirements, declined. 
Once all the questions and requirements for paper have been satisfied, 


the paper sent the New York Office with the appropriate recom- 
mendation. 


The work outlined above all done voluntarily and without remuneration. 
times the regular work those involved prevents immediate handling. 
may take much three weeks (more under special circumstances) for the 
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Chairman record receipt paper, choose reviewers and get copies 
the paper the reviewers. may require much three months for 
reviewer make his review and report the Chairman, depending his 
work load. (He asked within month, possible.) The review 
process can quite lengthy, particularly cases where only one copy the 
paper furnished. When extensive changes are required before publication 
can recommended, the time taken the author this work may vary 
months, longer, depending the time the author willing able 
spend the paper and the extent the revisions required. 

You, authors, reviewers, and prospective reviewers can speed and 
improve this process by, (a) preparing your paper meet all the re- 
quirements listed the inside front cover each issue the Journal; (b) 
signifying your willingness review papers your fields interest (see 
form back cover); (c) making prompt reviews papers sent you for 
review; and (d) making, expeditiously, changes comments clarify 
questions and requirements brought out the review. 


HYDRAULICS DIVISION ACTIVITIES DUE FOR FURTHER INCREASE 


The Division Newsletter October 1958 describes briefly the report made 
March the ASCE Task Committee Water Resources. recommend- 
the report, ASCE Water Resources Coordinating Committee was ap- 
pointed December 1958. Those appointed the Committee are: Finley 
Laverty, Chairman, Joseph Tiffany, representing the Hydraulics Division, 
William Donnan, the Irrigation and Drainage Division, George Vencill, 
the Power Division, Richard Kennedy, the Sanitary Engineering Division, 
and Roger Gilman, the Waterways and Harbors Division. Don Reynolds 
was assigned committee secretary. 

The purpose the Coordinating Committee coordinate and supervise 
the new organization and activities which are expected result from the 
very significant study made the Task Committee Water Resources. 
the next step its activities the Coordinating Committee has asked each 
participating Division undertake study its own organization, related 
the recommendations the Task Committee, and furnish recommen- 
dations the Coordinating Committee. Suggestions the Executive Com- 
mittee, Hydraulic Division, include request for specific proposals for ad- 
ditional committees and personnel man these committees. preliminary 
report from the Executive Committee was due early March 1959. 

The actions proposed the report the Task Committee indicate that 
implemented the Hydraulics Division the one most directly affected the 
field water resources. The new programs and studies that would become 
the responsibility the Hydraulic Division would extensive and im- 
portant approximately equal all other activities the Division. 


IMPROVEMENT PUBLICATIONS STANDARDS URGENTLY NEEDED 


The minutes various committee meetings well correspondence be- 
tween individuals dealing with ASCE publications contain many references 
the current lack proper recognition given the Division Journals ASCE 
prospective authors, readers, and library officials. The Journals their 
present form are looked upon interim measure pending publication 
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some the papers ASCE Transactions. Thus, the general membership 
deprived complete coverage some the best papers direct interest 

them. Authors, having assurance that their work will appear later 
ASCE Transactions, naturally give preference journals other recognized 
publications other professional societies. this manner some outstanding 
papers written ASCE members are lost insofar the records and the 
history the Society are concerned. Moreover, this lack recognition 
the Journals not only affects the attitude the author toward producing 
worthwhile paper, but also the attitudes those individuals who are asked 
review the paper well those who participate discussions the techni- 
cal material. 

appears therefore that action improve the standards publication for 
the Journals urgently needed. this end suggestions have been made 
the Publications Committee, ASCE, that early attention given converting 
the Journals into suitable and attractive media for publication professional 
papers, including the adoption publication standards that will recognized 
the Library Congress and other scientific libraries being appropriate 
for permanent inclusion their collections. 


INTERESTING PAPER HYDROLOGY PUBLISHED JOURNAL 
THE SANITARY ENGINEERING DIVISION 


The paper “Hydrology Urban Runoff” Tholin and Keifer 
being published the March Journal the Sanitary Engineering Division 
Proceedings Paper No. 1984. Members who are interested obtaining copies 
this paper can order using the coupon the “Proceedings Available” 
section the March issue CIVIL ENGINEERING. 


FOR YOUR CALENDAR 
ASCE Meetings 


May 4-8, 1959 
July 1-3, 1959 


October 19-23, 1959 
March 7-11, 1960 
June 19-23, 1960 
October 9-13, 1960 
April 10-15, 1961 
October 16-20, 1961 
February 1962 

May 1962 

October 15-19, 1962 


April 13-15, 1959 


June 15-19, 1959 


ASCE, Cleveland Convention 

Hydraulics Division, Fort Collins 
Conference 

ASCE, Washington, C., Convention 

ASCE, New Orleans Convention 

ASCE, Reno Convention 

ASCE, Boston Convention 

ASCE, Phoenix Convention 

ASCE, New York Convention 

ASCE, Houston Convention 

ASCE, Omaha Convention 

ASCE, Detroit Convention 


Non-ASCE Meetings 


ASME Hydraulic Conference, University 


Michigan 
American Society for Engineering 


Education, Pittsburgh, Pennsylvania 
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June 17-20, 1959 


June 18-20, 1959 


August 24-29, 1959 


September 9-11, 1959 


April, 1959 


National Society Professional Engi- 
neers, New York City 

ASME Applied Mechanics Conference 
Virginia Polytechnical Institute, 
Blacksburg 

International Association for Hydraulic 
Research, Montreal, Canada 

Midwestern Conference Fluid and 
Solid Mechanics, University 
Texas, Austin 


Deadline dates for Newsletter contributions: June issue April 15. 


August issue June 15. 


USE THE HYDRAULICS DIVISION NEWSLETTER 


You are urged continue use the Division Newsletter for announce- 
ments, inquiries, personnel news, committee reports, surveys and other items 
interest Division members. short note covering the highlights com- 
mittee meetings (including task forces) particularly requested. 


GUY ARBUTHNOT, Jr. 
Box 631 

Vicksburg, Mississippi 
Newsletter Editor 


P.S. Please fill out the form the next page now and return Mr. Ball. 
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REVIEWERS NEEDED 


Your help earnestly solicited obtaining willing and qualified members 
who are experts their fields review technical papers. The form below 
furnished for your convenience nominating yourself, and/or others for 
this important work. 
Reviewers are needed the following fields: 
Research 15. Tides and Waves 
Fluid Mechanics 16. Beaches and Estuaries 
Hydraulic Models 17. Ice and Permafrost 
Open Channel Flow 18. Water Measurement 
Closed Conduit Flow 19. Design 
Turbulence 20. Spillways 
Density Currents 21. Outlet Works 
Cavitation 22. Canal Structures 
Sedimentation 23. Gates and Valves 
10. Hydrology 24. Operation and Maintenance 
11. Floods 25. Analogs 
12. Ground Water 26. Instrumentation 
13. Water Supply 27. Other 


14. Meteorology 


Your name and address: Indicate fields interest number, 
thusly: (1), (2), 


| 
Name and address recommended Indicate Reviewer’s fields interest 


reviewers: number: 

Mail to: 


James Ball, Chairman 
Committee Publications 
Hydraulics Division, ASCE 

Bureau Reclamation 

Building 56, Denver Federal Center 
Denver 25, Colorado 
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